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ABSTRACT 

We demonstrate a novel technology that combines the power of the multi-object spec- 
trograph with the spatial multiplex advantage of an integral field spectrograph (IFS). 
The Sydney-AAO Multi-object IFS (SAMI) is a prototype wide-field system at the 
Anglo- Australian Telescope (AAT) that allows 13 imaging fibre bundles ("hexabun- 
dles") to be deployed over a 1-degree diameter field of view. Each hexabundle com- 
prises 61 lightly-fused multimode fibres with reduced cladding and yields a 75 percent 
filling factor. Each fibre core diameter subtends 1.6 arcseconds on the sky and each 
hexabundle has a field of view of 15 arcseconds diameter. The fibres are fed to the fiex- 
ible AAOmega double-beam spectrograph, which can be used at a range of spectral 
resolutions (i? = A/5A « 1700-13000) over the optical spectrum (3700-9500 A). We 
present the first spectroscopic results obtained with SAMI for a sample of galaxies at 
z « 0.05. We discuss the prospects of implementing hexabundles at a much higher mul- 
tiplex over wider fields of view in order to carry out spatially-resolved spectroscopic 
surveys of 10* — 10^ galaxies. 

Key words: instrumentation: spectrographs - techniques: imaging spectroscopy - 
surveys - galaxies: general - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Galaxies are intrinsically complex with multiple components 
and varied formation histories. This complexity is the pri- 
mary reason that unravelling the physics of galaxy forma- 
tion and evolution is so challenging. Galaxies are made up 
of baryons confined to dark matter haloes, and often have 
multiple distinct kinematic components (e.g. bulge and/or 
disc). There are complex interactions between the stars, gas, 
dust, dark matter and super-massive black holes. These can 
lead to both positive and negative feedback on the formation 
rate of stars. 
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Experimental efforts to address galaxy formation have 
generally taken two directions. First, galaxy imaging and 
spectroscopic surveys have progressively moved to higher 
redshift, in an attempt to directly observe galaxy evolution 
and formation. This approach has had much success, plac- 
ing strong constrai nts on the evo lution of the global star 
formation rate (e.g. lHopkins fc Beacom 200 6). unveiling the 
strong evo lution of black hole accret i on over most of cosmi c 
time (e.g. ICroom et all |2004 l2009l : [Richards etlll bood ). 
tracing the evolutio n of galaxy size and morphology (e.g. 
iDressler et aTlllQQTl ) and much more. 

The second approach has been to expand our view in 
wavelength rather than cosmic time. The physical processes 
occurring in galaxies cause emission over the entire range of 
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the electromagnetic spectrum. In order to have a full pic- 
ture of galaxy formation, a multi-wavelength approach is 
vital. This has been made more achievable with recent gen- 
erations of satellites covering the X-ray, ultra-violet, mid- 
and far-infrared. While the spectral energy distributions of 
stars tend to peak in the optical or near-infrared, obscu- 
ration and reprocessing by dust generates strong mid- and 
far-infrared emission. Young stars (when not obscured) are 
most directly traced in the ultraviolet, while black hole ac- 
cretion can generate radiation from the radio to X-ray and 
gamma-ray bands. Surveys such as the Gala xy And Mass 
Assembly (GAMA) survey llDriver et al.ll201lh and the Cos- 
mic Evolution Survey rCOSMOS: IScoville et al1l2007l ) show 
the value of this multi- wavelength approach. 

The third route, and the one that we address in this 
paper, is to focus on spatially resolving galaxies; in particu- 
lar, obtaining spatially resolved spectroscopy. Optical spec- 
troscopy allows us to measure a wide range of parameters 
including current star formation rates (e.g. via Ha), gas 
phase metallicities, stellar ages, stellar metallicities, black 
hole accretion rates, ionization structure and extinction due 
to dust (e.g. via the Balmer decrement). The major spcctro- 
scopi c surveys to date have u sed a single fibr e (CoUcss et aL i 



2001 ; lAbazaiian et al.ll2003h or single slit (jLe Fevre et al.l 



20051 : iDavis et al.ll2007f ) on each object, and so obtain just 



one measurement of these parameters for each galaxy. More- 
over, these measurements may not be representative of the 
galaxy as a whole, but biased, depending on where the aper- 
ture is placed. This fundamental problem is addressed by in- 
tegral field spectrog raphs (IFS). In the last decade, projects 
such as SAURON H aeon et al.l I2001I ) have demonstrated 
the power of integral field spectroscopy to capture a range 
of key observables that are simply not available to single- 
aperture surveys. As well as studying the properties listed 
above in a spatially-resolved context, obtaining gas and stel- 
lar kinematics over an entire galaxy enables us to separate 
dynamical components, measure dynamical mass, examine 
the impact of winds and outfiows, and discover merging sys- 
tems via dynamical disturbances. 

Integral field spectroscopy has almost exclusively been 
limited to single-object instruments, meaning that it 
is time-consuming to build large samples. The largest 
current data set, using the SAURON system on the 
William He rschel Telescope, con tain ~ 260 objects 
fATLAS-3D:ICappellari et al.ll2011al '). The CALIFA project 
jSanchez et al.ll201ll ) aims to target 600 objects with the 
PMAS integral field unit (IFU) on the Calar Alto Telescope 
using ~200 nights of telescope time. The only multi-object 
integral field spectrograph cu rrently available is that on 
the VLT FLAMES instrument l|Pasquini et al.ll2002l '). which 
contains 15 IFUs, each with 20 spatial elements of size 0.52 
arcsec and a total field of view (FoV) of 2 x 3 arcsec. This has 
enabled m easurements of the TuUy-Fisher relation at ~0.6 
l|Yang et al. 2008) as well as a number of other projects. 
However the small FoV of each IFU, combined with the high 
spectral resolution (>9000) and associated narrow wave- 
length range limits its applicability for large-scale surveys. 

Astrophotonic technology ( Bland-Hawthorn fc Kerr] 

I2OO9I ) is now opening the way to new instrumentation 
that can address the need for hi ghly multiplexed inte- 
gral field spectroscopy. H exabundles (|Bland-Hawthorn et al.l 
l201ll : [Brvant et aLll201ll 'l are optical fibre bundles where the 



cladding has been stripped from each fibre to a minimum 
over a short length (~30mm) and the fibres then gently 
fused together at the input end to provide an IFU ('^Imm 
aperture) with high filling factor. These can then be used in 
conventional multi-fibre spectrographs. 

In this paper, we report on the Sydney-AAO Multi- 
object Integral field spectrograph (SAMI), the first fully 
operational demonstrator instrument to use hexabundles. 
SAMI has 13 IFUs that can be positioned anywhere over a 
1 degree diameter field of view. In Section [21 we discuss in 
detail the scientific rationale for such an instrument, along 
with some practical considerations regarding sensitivity. In 
Section [31 we describe the SAMI instrument in detail. In 
Section[31 we outline the observations carried out during the 
commissioning of the instrument, the results from which are 
discussed in Section [SI Our conclusions, and goals for the 
future, are laid out in Section [SI 



2 SCIENTIFIC RATIONALE 

In this section we outline the key scientific drivers for an 
instrument such as SAMI. The fundamental question at the 
heart of this work is: how did the galaxy population we see 
around us today come about? This requires us to under- 
stand the physical processes that occur as galaxies form and 
evolve. The galaxy population we see today has some very 
distinctive features that need to be explained. 

One of the most fundamental is the separation of galax- 
ies into a bimodal distribution accor ding to colour (e.g. 
IStrateva et al.l200ll : [Baldrv et al1l2004 ). A galaxy's colour is 
primarily related to the presence of ongoing star formation. 
The second key feature differentiating galaxies is morphol- 
ogy. There is a strong correlation between colour and mor- 
phology, with galaxies lying along the 'red sequence' being 
mostly passive systems with elliptical/spheroidal morphol- 
ogy, while galaxies inhabiting the 'blue cloud' are mostly 
dominat ed by discs, although this is not always stri ctly the 
case (e.g. [Masters et al.ll201(]| : ISchawinski et al.ll20 09!') . While 
they are related, there is not a strict one-to-one relationship 
between morphology and colour. A clearer understanding 
of galaxies can be obtained if they are considered as be- 
ing made up of distinct morphological components (discs, 
bulges and pseudo-bulges) that result from different for- 
matio n processes and evolutionary histories (jPriver et al.l 
l2007h . The intrinsic properties of these structural compo- 
nents are more uniform than those of the galajcies they 
compose. Their formation pathways are also quite differ- 
ent, with true bulges built up by violent mergers, discs 
from gas accretion, and pseudo -bulges from secular evolution 
(jKormendv fc Kennicuttll2004l ). Disentangling these various 
modes is complex, but can be materially aided by the fact 
that the structures formed have different kinematic proper- 
ties as well as different star formation histories. 

In order to ascertain the physics involved, we need to 
determine the answers to a number of questions that broadly 
fall into four categories: (i) how does galaxy mass and angu- 
lar momentum build up? (ii) when, where and why does star 
formation occur? (iii) when, where and why does black hole 
accretion occur? (iv) how are galaxies fuelled and what is 
the role of feedback? We will discuss each of these issues in 
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turn, although there is of course significant overlap between 
them. 



2.1 The build up of mass and angular momentum 

The standard picture of galaxy formation has gas cooling to 
form a rotationally supported disc within a cold dark mat- 
ter halo (White & Rces 1978). While this picture is broadly 
accepted, feedback and interactions provide major compli- 
cations which are not yet fully understood. 

The scaling relation between circular v elocity and stel- 
lar in ass (The TuUy-Fisher Relation, TFR; ITuIIv fc Fisheil 
I1977I ) for disc galaxies provides a tight constraint for galaxy 
formation models. The circular velocity depends on the ra- 
tio of disc mass to halo mass, the dark matter halo profile 
and the dimensionless spin parameter, A dPe cblcs 1963). The 
largest TFR samples (e.g. ISpringob et al.l 120071 ') currently 
contain ~ 5000 galaxies with either long-slit spectroscopy 
or spatially unresolved H I velocities. If IFU observations 
can be made out to large enough radii (typically ~ 2.2 disc 
scale lengths), then they provide substantial advantages in 
allowing a clearer picture of distorted kinematics and incli- 
nation. Circular velocities can be compared to the results of 
galaxy lensing to constrain the dark matter halo profile and 
examine evi dence of contraction of the halo in res ponse to 
the baryons (iDutton et al.ll2010l : iReves et al.ll201lf ). 

The stellar and emission line kinematic data that SAMI 
can provide will allow dynamical mass estimates within max- 
imum radius pr obed, using techn iques such as anisotropic 
Jeans modelling l|Cappellarill2008l ) . In general it is not possi- 
ble to determine total mass because of the uncertainty of the 
dark matter halo parameters. Even for the Milky Way, where 
many halo stars can be used to probe the outer halo, the to- 
tal ga laxy meiss is uncertain to a factor of two (jSmith et al.l 
[2003) • Detailed dynamical techniques are in stark contrast 
to estimates of dynamical mass which simply take a single 
velocity dispersio n of the gala^xy (e .g. from a single fibre 
observation), (e.g. iTavlor et aLll2010l ). 

It is relatively easy to extract a rotation curve v{r) 
from t he observed data if the kine matics are fairly well or- 
dered (|Stavelev-Smith et all [19901 ). This can then be used 
to provide dynamical information about the galaxy, partic- 
ularly if baryonic information is brought to bear. However 
if the kinematic axes are misaligned with the photomet- 
ric axes, this is often a signature of streaming motions due 
to a bar or an oval distortion. A dynamical mass can still 
be derived, although the increa sed number of free parame- 
ters makes this more uncertain (IStavelev-Smith et al.|[l990l : 
iQuillen. Frogel fc Gonz alez 1994). Any deviations from ro- 
tational symmetry are important in their own right. It is 
often very difficult to see the presence of bars, particularly 
in highly-inclined disc systems. But bars are often betrayed 
by the inner twists of the isovelocity contours. Warps are 
more easily de tected on lar ge scales and almost always in 
HI kinematics l|Briggslll990l ): however, the same effects can 
now be seen in dee p observations of the diffuse ionized gas 
in th e outer disc l|Christlein. Zaritskv fc Bland-HawthornI 
l2010h . The physical cause of kinematic distortions can be 
examined by large IFU surveys which probe a variety of 
galaxy parameters. For example, are distortions more likely 
in high density regions due to dynamical interaction. 

Integral field spectroscopy also enables studies of stellar 



kinematics that describe the observed projected stellar an- 
gular momentum per unit mass of galax i es, not possible with 
other techniques ijEmsellem et al.ll2007l . I2OIII : iBrough et al.1 
l201lh . This measurement enables the separation of early- 
type galaxies into fast and slow rotators which are thought 
to have very different formation paths. A great success of 
the SAURON and ATL AS-3D projects (|Bacon et al.ll200ll : 
ICappellari et al.]|2011ah has been the discovery that most 
early type galaxies have significant rotation, with only ~ 14 
percent (predom i nantly at high mass) being slow rotators. 
ICappeUari erahl (|2011bh have used ATLAS-3D to demon- 
strate a kinematic morphology- density relation, which shows 
a smooth transition of spirals to early type fast rotators with 
increased density, and massive slow rotators only inhabiting 
the highest density regions. SAMI would allow such studies 
to be expanded to probe greater dynamic range in environ- 
ment and examine a such relations as a function of stellar 
mass. 

The rate of dark matter halo merging can be ac curately 
estimated from simulations (|Fakhouri fc Mai l2008t ) . Kine- 
matic information from integral field spectroscopy can be 
used to differentiate betw een quiescent galax ies and those 
undergoing a merger (e.g iShapiro et al.ll2008l) usmg proce- 
dures such as kinemetry (iKrainovic et al.l 20061 ). Until now, 
this type of analysis has largely been limited to high red- 
shift (where the merger rate is expected to be higher), but 
with samples of 10'' or more galaxies, statistically mean- 
ingful estirMtes_ofmerge£rates can be made at low red- 
shift. iFakhouri fc Mai l|2008l ) predict that the halo merger 
rate at z ~ should be ~ 0.05 halo"'^ Gyr "'^ for major 
mergers (with mass ratios < 3/1). Integral field observa- 
tions provide a complementary approach to studies which 
focus on the number of close pairs to estimate merger rates 
(e.g. iDe Propris et al.ll2007^ , as they probe very different 
phases of the merger process. It is also be possible to look 
for weaker dynamical distur bances in discs due to repeated 
minor mergers/interactions jZaritskv fc Rixlll997l ). In this 
case, it is very advantageous to extend the study to large 
radius as the effect of tidal perturbations scales as ^ r"^. In 
this case the challenge is to achieve sufficient sensitivity at 
large radius (e.g. >1.5 — 2 scale lengths). 

Spin angular momentum from galaxy kinematics can 
directly probe the formation of the large scale structure of 
the universe and galaxy formation. The spin of galaxy discs 
provides an approximation of the spin of the galaxy's dark 
matter halo (|Sharma fc Steinmetzi2005 ) , which is coupled to 
the large scale structure. Early in a dark matter halo's life, it 
experiences torques from the surrounding density landscape. 
The spins of dark matter haloes today retain some memory 
of that landscape, so spin is intrinsically linked with the 
large scale structure. This link may be examined in N-body 
simulations and observations by measuring the distribution 
of inferred dark matter halo spin magnitude (|Berta et al.l 
2008) and by the relative orientation of galaxy spin direc- 
tions with each other and with the large scale structure. 

N-body simulations do not predict a strong align- 
ment between the spins of neighbouring haloes, although 
an alignment between galaxies has bee n detected in the 
TuUy catalogue of 12,122 nearby spirals (|Pen. Lee fc Selia^ 
I2OOOI ). Simulations and theory predict an alignment be - 
tween halo spin and the tidal field l|Lee fc Peiij l2000l ). 
and an alignment with features in the tidal field like fil- 
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aments (IZhang et al l 120091 '). sheets l|Leel |2004 ) and voids 
l|Brunino et al.l2007 ). Any kind of alignment is predicted to 
be very weak, however, so could only be seen in large scale 
galaxy surveys. There have been detections of spin align- 
ments in the large scale structure reconstructed from imag- 
ing su rveys (|Lee fc Erdogdull2007l : IPaz. Stasvszvn fc Padillal 
I2OO8I ). using inferred galaxy spins from disc shapes. A 
tentative detection of spin alignment with filaments was 
found using the inferred spin of only 2 01 galaxies around 
voids jTruiillo. Carretero fc Patirill2006D and 70 galaxies in 



filaments ( Jones, van de Weveaert fc Aragon-Calvd l2010h . 
picked from the large scale structure of SDSS. Discrepancies 
between the results found from dark matter simulations and 
observations indicate a difference in the way that galaxies 
and dark matter haloes obtain and keep their spin. A large 
survey of direct spin measurements could reveal whether 
galaxies exhibit the same spin behaviour as dark matter 
haloes, and show how galaxy spin is linked to the large scale 
structure. 



2.2 When, where and why does star formation 
occur? 

Much recent observational and theoretical work has focussed 
on how blue galaxies can have their star formation quenched, 
moving them onto the red sequence. Red sequence galax- 
ies are preferential ly found in denser environments (e.g. 
iBlanton et"aLll2005l ). and star formation is also c learly sup- 
pressed at high density (e.g. iLewis et al.ll200^ . This im- 
mediately suggests environmental factors play an important 
role. When a galaxy falls into a clust er, the ram pressure 
from the dense intergalactic medium (|Gunn fc GottlllOT^ ) 
may expel the gas from the disc, removing the fuel re- 
quired for further star formation. There are several observed 
examples of this in r i ch cl usters (e.g. iRandall et al.l [2OO8I : 
ISun. Donahue fc Voitl|200^ . In moderately dense regions, 
such as galaxy groups, it may be that ram pressure will 
leave the disc intact, but can still remove gas from the halo 
of the galaxy. The halo provides a reservoir of gas which can 
replenish the disc. 

Without the halo gas, the star formation will decline 
and then cease as the disc gas is consumed, leading to a tran- 
sition to the red sequence, in a process known as strangu- 
lation l|Larson. Tinslev fc Caldwell|[l980l ). Simulations sug- 
gest that this process can be efficient at re moving halo 
gas, even in small a nd /or compact groups (|Bekkil [2OO9I : 
iMcCarthv et al.l2008l '). but there is little direct experimental 
evidence that this is the case. Indirect evidence does point 
to some pre- processing of galaxies in groups before they fall 
into clusters (jBalogh fc McGe3l2O10l ). but the physical pro- 
cess driving this is not clear. Direct galaxy-galaxy interac- 
tions are also expected to play a critical role , with major 
galax y mergers triggering star formation (e.g. ^Elli son et al.l 
|2008| ) and transforming the morphology of galaxies, al- 
though the fraction of galaxies undergoing major merg- 
ers (i.e. those with mass ratios of 3:1 or less) i n the lo- 
cal Universe is smaU (e.g. iPatton fc Atfieldlbooi ). On the 
other hand, dwarf star-forming galaxies in the local Uni- 
verse are often found interacting with low-luminosity ob- 
jects or diffuse H I clouds,. This appears to be the trig- 
gering mechanism of their intense sta r-formation activity 
l|L6pez-Sanchez fc Estebaiil I2OO8I . |2009| ). although only de- 



tailed mul ti-wavelength observations are able to reveal thes e 
processes (|L6pez-Sanchej|2010l : iLopez-Sanchez et al.ll201ll ) 
However, environment is only one factor. Feedback from 
star formation and accretion onto super-massive black holes 
provides an internal mechanism for transformation. This 
feedback provides a solution to the mismatch of the theoret- 
ical dark matter halo mass function and the observed stel - 
lar mass function (e.g. iBaldrv. Glazebrook fc Driver! [2OO8I ) 
by heating and/or expelling gas in both low mass (via 
star for mation) and high mass (via black hole accretion) 
haloes (|Cattaneo et al. I I2OO6I : lBaldr~ Glazebrook fc Driver] 
[200^. Extreme outbursts of star formation or black hole 
accretion may be tr iggered by mergers or interactions (e.g. 
iHopk ins et al II2O O8 '), making a link between internal and en- 
vironmental effects. Once the burst is over, another mecha- 
nism is needed to suppress continued star formation. The 
best suggestion for this is mechanical feedb ack from jets 
emitted by super-massive black holes (e.g. ICroton et al.l 
I2OO6I ). but this only appears to be efficient in massive galax- 
ies. 

As well as these active processes, the environment has 
an indirect influence via formation age. Galaxies in high den- 
sity regions form earlier and so have had more time to evolv e 
(e.g. lKaiseij|l984 iBardeen et al. l l 19861 : iThomas et aLll2005l ) . 
In the absence of other effects, we would then expect to see 
galax;ies in high density regions having older stellar popula- 
tions. 

Disentangling these varied influences on galaxy forma- 
tion is far from trivial. However, studies of the spatial dis- 
tribution of instantaneous star formation rates, integrated 
star formation (via stellar population ages) and metallicity 
(both gas and stellar) provide considerable insight. Impor- 
tantly, ram-pressure removal of gas implies that the trun- 
cation of star formation is an outside-in process (e.g. lBekkil 
l2009l : lKapferer" et ahllioogl ). Gas is preferentially removed in 
the outer parts of galaxies, which are less gravitationally 
bound. This may be a short-lived feature of the galaxies 
in dense environments if the gas is eventually completely 
removed. Alternatively, stripping can occur over the life- 
time of a galaxy if the gas is puffed up by the internal 
star formation; in this case, eve n a rarifled external medium 
can r emove gas from the galaxy l|Nichols fc Bland-HawthornI 
I2OIII ). Globally, the expectation would be that galax- 
ies form inside-out, and this implies age and metallic- 
ity gradients which are obser v ed (e.g. Sha ver et al. 'l983|; 
Vila-Costas fc Edmundi Il992l: [ Steinmetz fc Muellen ,1991 



Chiappini. Matteucci fc Grattonlll997f r 

One approach that has been explored in some detail as 
an alternative to spatiall y-resolved spectroscopy is the so- 
called 'pixe l-z' technique (|Conti et al.ll200^ : IWelikala et al.l 
|2009| . |2008| ). This approach, analo gous to photometric red- 
shifts ('photo-z'), uses a library of template spectral energy 
distributions (SEDs) to flt the observed optical and near- 
infrared colours of individual pixels within resolved galaxy 
images. This technique has been used with some success to 
explore t he environmental dep endence of star formation in 
galaxies. IWelikala et al.l (|2008l ) find that, globally, the sup- 
pression of star formation in high density environments (e.g., 
Lewis et al., ,2002 : .Gomez et al. 2003) seems to occur pri- 
marily in the most strongly star-forming population, and 
to be ev idenced by a s uppression in the inner regions of 
galaxies. IWelikala et al.l (2009) demonstrate that this effect 
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seems to hold independently for both early- and late-type 
galaxy populations, and that the suppression in star forma- 
tion cannot be explained solely by the w ell-known density- 
morphology relation (e.g.. lDressleij[l980l ). There are signif- 
icant limitations, however, to the 'pixel-z' approach. These 
are related to implicit assumptions made by the technique, 
such as each pixel being represented by an isolated single 
stellar p opulation with a simple exponential s tar-formation 
history (jConti et al.l |2003': 'Welik ala et al.l bOOSY As a result, 
the method cannot effectively measure the instantaneous 
star formation rate, which can be traced spectroscopically 
(e.g. by Hof emission). 

IFU spectroscopy allows us access to both current star 
formation (via emission lines) and integrated star formation 
history (via stellar age and metallicity). Examining the ra- 
dial dependence of the mean stellar age and metallicity gra- 
dients tells us when and where the stars formed in these 
galaxies, along with a fossil record of the galaxy merger 
history. The mean stellar age is effectively a luminosity- 
weighted integral of the star formation history whilst the 
stellar metallicity g radient provides an indication of its 
merging history fe.g. iKobavashil [20041 ). If environmental ef- 
fects are responsible for the cessation of star formation, then 
we would expect red sequence galaxies to have younger cen- 
tral ages with past major mergers sign-posted by shallow 
negative metallicity gradients (i.e. lower metallicities in the 
outskirts; lKobavashil2003 : lBrough" et al .l2007l : ISpolaor et all 
I2OO9I ). 



2.3 When, where and why does black hole 
accretion occur? 

A full picture of the physical processes involved in the fu- 
elling of accretion onto super-massive black holes, resulting 
in the phenomenon of an active galactic nucleus (AGN), 
is still elusive. It is now known that most galaxies contain 
super-massive black holes, with typic al masses a million to 
a billion times that of th e Sun (e.g. iGebhardt et al.l I2OO0I : 
iFerrarese fc Merritt|[200ol ') . The mass of the black hole corre- 
lates well with the mass ( or velocity disper sion) of the bulge 
or spheroidal component l|Tremaine et al.|[2002 l in a galaxy. 
This implies an intimate connection between the build up of 
stellar mass in galaxies and their super-massive black holes. 

The nature of the connection betwee n star forma- 
tion and AGN has long been debated (e.g. ISanders et al.l 
Il988t ). although a resolution remains elusive. The most 
luminous AGN (i.e. bright quasars) require >10'' solar 
masses of gas deposited in th eir central regions on time- 
scales of ~10^-10* years (e.g. ICroorn et^a l. 2005 ), requir- 
ing m ajor galaxy- wide perturbations ( Hopkins fc Hernguistl 
[2009|), and no doubt triggering substantial star formation. 
In contrast, low- luminosity AGN require relatively small 
amounts of gas, which can be supplied by internal stochas- 
tic processes, such as the ac cretion of cold molecular clouds 
l|Hopkins fc Hernauistll2009l ). 

The AGN population evolves strongly with cosmic 
time. This is most clearly seen in evolution of luminous 
quasars from z = to z ~ 6 (|Richards et al.l l2006l : 
ICroom et al.|[2009l ). which have a strong peak in space den- 
sity at 2 ~ 2 — 3. Lower luminosity AGN show less severe 
evolution, and peak in space density at lower redshif t (e.g. 
iHasinger. Mivaii fc Schmidtll2005l : ICroom et aLll2009l ). This 



downsizing appears q ualitatively similar to that seen in the 
formation of galaxies l|Cowie et al.l 1 199^ . In the local Uni- 
verse, accretion onto black holes is dominated by systems 
with low black hole masses, < 10* M0, but which are typi- 
cally accre ting within an order o f magnitude of the Edding- 
ton limit (|Heckman et al.ll2004l ). In contrast, higher mass 
black holes in the local Universe are accreting with a char- 
acteristic timescale substantially longer than a Hubble time. 

Luminous AGN (as measured by their [O III] A5007 
luminosity) are found to have younger stellar populations 
(jKauffmann et al1l2003l ). as measured by their Dn(4000) and 
H5 line indices from SDSS spectra. It is not clear whether 
this younger stellar population is centrally concentrated, or 
is distributed more widely across the host galaxy, as the 
SDSS fibres subtend a physical scale of ~ 5 kpc a.t z — 0.1. 
IFS data has the ability to explicitly examine the distri- 
bution of star formation and stellar population ages across 
galaxies and investigate whether these are related to accre- 
tion rate. 

Lower luminosity AGN tend to have spectra typical 
of Low-Ionizati on Nuclear Emission-line Regions (LINERS; 
iHeckmanirigSd ). While there is evidence that these form 
a continuous sequence with Seyferts (|Kewlev et al.l [20061 ). 
there substantial evidence that LINER like emission is ex- 
tended and may be p owered by ionization from AGB stars 
l|Yan fc Blantonll201ll ). A large galaxy survey with spatially 
resolved spectroscopy could resolve this issue. If LINER 
emission is not due to a central AGN in most cases this 
would require substantial re-interpretation of recent work 
on low redshift AGN. At an even more fundamental level, 
the fraction of galaxies whic h host AG N is only well de- 
fined locally (i.e. 10s of Mpc; [Hc1I20o1 ). Approximately 40 
per cent of these very local galaxies host AGN. At larger 
distances contamination from off nuclear emission increas- 
ingly reduces the sensitivity to weak nuclear emission lines. 
A large IFS survey would enable apertures of fixed metric 
size to be defined, enabling robust AGN rates as a function 
of redshift to be determined. 

While there is a good theoretical basis for galaxy merg- 
ers t riggering luminous AGN (e.g. iHopkins fc Hernguistl 
|2009| ). the evidence for this is mixed. An alternative path- 
way is via violent disc instability in self grav itating discs 
(e.g. iBower et al. I I2OO6I ; iBournaud et al.] 120111 ). Disc insta- 
bility is at least likely to play a role at high redshift where 
cold streams can dominate the mass accretion onto discs 
l|Dekel et al.ll2009l ). Indeed, a large fraction of high-redshift 
star f orming galaxies appea r to be discs (e.g. iGenzel et al.l 
l2008l : IWisnioski et al.ll201ll ). rather than mergers. In low 
redshift samples AGN activi ty is not enhanc ed by the pres- 
ence of a n earby companio n ijLi et al.ll2008lj ) , while star for- 
mation is l|Li et al.|[2008al ). This is a somewhat surprising 
result, and may point to a difference in timescale between 
the onset of star formation and the AGN, with the AGN oc- 
curring later, after the merger has taken place. Indeed, IFS 
observations of local galaxies with observed outfiows demon- 
strates that the AGN ti mescale is significantly longer tha n 
the starburst timescale (|Sharp fc Bland-HawthornI I2OIOI ). 
Large scale IFS observations can directly address the issue 
of AGN fuelling by examining the kinematic properties of 
AGN hosts, and searching for evidence of disc instability 
and/or merging. In this regard it will be particularly impor- 
tant to span a range of accretion luminosities in order to 
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examine whether there is a change from secular evolution at 
low accretion rate to mergers at high accretion rate. 

2.4 Feeding and feedback 

The accretion of gas onto galaxies remains a largely un- 
solved problem. Whether gas enters the galaxy potential in a 
hot p hase and then cool s down ([ Binncv, Nipoti &_Pratcrnali 
I2OO9I ). as a warm rain l|Bland-Hawthorn et allboOTi ). or as 
an H I complex like the h igh velocity clouds in the Galactic 
halo H ancisi et al.ll2008h . or all of th e abov e, is unclear. Af- 
ter a review of the evidence, iBinnevI l| 19921 ') concluded that 
the outer warps of HI discs were some of the best evidence of 
ongoing disc accretion. Once the gas settles into the galaxy, 
the outstanding issues are how the gas feeds into the nuclear 
regions and, in particular, onto a central black hole. 

In a survey of 10^ - 10'' galaxies, SAMI offers 
the prospect of studying nuclear activity (AGN, star- 
burst, LINER) and star formation within the context of 
the extended galaxy. The kinematic signatures of outer 
warps and i nner bar streaming are rel atively easy to pick 
out in HI llStavelev-Smith et al.l Il990l) or in ionized gas 
l|Christlein. Zaritskv &: Bland-HawthornI bOloD . Thus, we 
can now directly associate this activity with large-scale disc 
disturbances, assuming these exist. Traditionally, in large 
galaxy surveys, the association of activity and dynamical 
disturbances is made from the proximity of galaxies in po- 
sition and redshift sp ace (e.g. iNikolic. CuUen fc Alexander! 
I2OO4I : iLi et ahlliboSah . With fuU kinematic information a 
much more direct determination of the triggering of activity 
will be possible. 

An alternative approach is to study the impact o f the in- 
ner d i sc on the extended pr operties of galaxies (e.g. iMartinI 
1 19981 : iKing fc Pounds! [20o3 ) . Jets carry energy, and winds 
carry gas and metals, far from the nucleus. In a recent 
integral field study of ten AG N and starburst galaxies, 
ISharp fc Bland-HawthornI (|2010l ) find that starburst winds 
are largely shock ionized, while AGN winds show the hall- 
mark of photo-ionization by the accretion disc, clearly indi- 
cating that the starburst phenomenon is very short. For one 
of the objects observed on the SAMI commissioning run, we 
see the ionization characteristics typical of nuclear activity 
for gas off the plane of an inclined disc (see Section [5.41 and 
Fogarty et al., in prep.). A large-scale wind is confirmed 
by the broad emission-line profiles along the minor axis. 
This is a remarkable testament to the power of spatially- 
resolved kinematic and ionization information. In the full 
survey, SAMI is likely to uncover hundreds of new outfiow 
sources connected either to nuclear activity or inner disc star 
formation. There is an even rarer class of galaxie s with disc- 
wide winds that SAMI will also inevitably add to (|Stricklandl 
I2OO7I ). 

A number of edge-on spiral disc galaxies hav e ver- 
tically extende d ionized gas in their haloes ( Rai3 1 19961 : 
iDettmaij 1 19921 '). The Reynolds Layer in our Galaxy, re- 
cently mapped by the WHAM Ha survey telescope 
l|Madsen. Haffner fc Reynolds! 2006!). is a good example of 
this phenomenon (see also 'Gaenslor et al.! l2008l ). The ion- 
ization characteristics of this gas does not appear to be con- 
sistent with any known mechanism (i.e. UV photo-ionization 
by hot young stars, radiation from old supernova bubbles, 
shocks from supernovae, cosmic ray heating, or radiation 



pressure on dust grains in the disc). The gas may arise 
from some kind of disc-wide interaction between the disc 
and the hot halo, presumabl y driven b y processes related 
to star formation in the disc (|Coxll2005! l But there is also 
the prospect that some of this gas is related to warm gas 
accretion onto the disc (|Bland-Hawthorn!l2009l ') or involved 
in a large -scale circulation or recycling of gas through the 
halo (Mar inacci et al.ll20lih . 

As with the study of galaxy winds, a large SAMI survey 
has the potential to greatly increase the sample of known 
galaxies with vertically-extended warm discs. With a larger 
sample, it will be possible to correlate the presence of these 
discs with the disc star formation rate, nuclear activity and 
galaxy mass. 



2.5 Limitations of current spectroscopic surveys 

Historically, telescopes were used to observe one source 
at a time. But with technical advances in optical fibres, 
it was realized in the early 1980s that many sources 
could be observed simultaneously across the telescope 
focal plane by precisely positioning fibres in the field 
(|Barden. Ramsev fc Tru"5 Il98ll : iGrav! Il983!) . This led to 
an explosion in wide-field spectroscopic surveys, notably in- 
cludi ng the 2-degree Fi eld Galaxy Redshift Survey (2dF- 
GRS:IColles8 et al.ll200ll ). 2-degree Field QSO Survey (2QZ 
Croom et al.l l2004 ). 6-degree Field Galaxy Survey (6dFGS 
Jones ct al. 2009) and Sloan Digital Sky Survey (SDSS 
York et al.i i2QQd ) amongst several others. Between them, 
such surveys have obtained spectra for approximately 1.5 
million extragalactic targe ts. New instrum ents recently com- 
missioned (e. g. LAMOST:[Su et al. l ll998! ) or in construction 
(e.g. VIRUS: lHill et al.ll2004 ) are able to observe thousands 
of sources at a time. 

Considerable advances have been made possible by the 
2dFGRS and SDSS, which use a single optical fibre per 
galaxy. However, with fibre diameters of 2 and 3 arcseconds 
respectively, these projects sample less than half the light 
from a galaxy at the median distance of the surveys. These 
single apertures limit the surveys in two ways. 

First, it is impossible for single-fibre surveys to mea- 
sure spatially varying spectral properties, which prohibits 
the study of crucial observables such as kinematic merger 
rates, galaxy rotation and dynamical mass, star formation 
gradients, metallicity gradients, age gradients and detection 
of galaxy winds and/or outfiows. 

The second limitation is that with single-fibre spec- 
troscopy the measured signal depends on many things: (i) in- 
trinsic properties, like source luminosity, size and distance; 
(ii) atmospheric conditions, particularly seeing; (iii) instru- 
mental properties, like fibre aperture size and positioning 
accuracy, and optical focus over the field; (iv) telescope prop- 
erties, such as pointing and guiding precision; and perhaps 
other effects. Many published papers make the mistake of as- 
suming the surveys provide spectrobolometry (i.e. the spec- 
trum of the total light output by the source) rather than 
the spectrum from an (often ill-defined) spatial sample of 
the source. The inherent dangers of aperture effects have 
long been known in astr onomy, but have often been under- 
appreciated or ignored. lEUis et aL ] l|2005!, see their Fig. 8) 
clearly demonstrate that the single- aperture fibre spectra 
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from a typical galaxy survey may be only weakly correlated 
with the photometric classification. 

Aperture biases can manifest themselves in two ways. 
The fibres subtend an increasing linear size with increasing 
distance of the galaxy, potentially causing spurious evolu- 
tionary effects (and spurious luminosity-dependent effects 
in a flux-limited sample). Second, important galaxy prop- 
erties, such as star formation rate and metallicity can have 
substantial gradients, meaning that observations of just the 
central regions are not representative of global values (e.g. 
iKewlev fc Ellisonll2008l '). 

2.6 Why multiplexed IFUs? 

Integral field spectroscopy allows us to gather data on key 
observables that are simply inaccessible to single apertures. 
IFUs are usually a single monolithic array of lenslets, with 
the light fed to a spectrograph via optical fibres (or alter- 
natively an optical image slicer). As a consequence, they 
can only target one object at a time and so IFU surveys 
have typical ly only targeted a few tens of galaxies in specific 
classes (e.g. lPracv et al.ll2009l : fEmsellem et al.ll2007l ). 

The process of galaxy formation and evolution is in- 
herently complex, with the observed galaxy properties de- 
pending on a large number of parameters, such as host halo 
mass, stellar mass, merger history etc. As well as the multi- 
dimensional nature of galaxy properties, there is inherent 
stochasticity in the process. This is at least in part due to 
our inability to accurately trace the formation history of in- 
dividual objects, but also derives from inherently non-linear 
physics such as that involved in the collapse of molecular 
clouds to form stars. 

The multi-dimensional nature of the galaxy population, 
combined with this stochasticity (which at some level could 
be considered as extra hidden parameters) , means that large 
surveys are required to extract the key relations between 
physical properties. The success of surveys such as 2dFGRS 
and SDSS has in large part been due to their ability to 'slice 
and dice' the galaxy distribution and still have statistically 
meaningful samples in each bin of the parameter space. 

The need to extend integral field observations to large 
samples is well understood in the above context, and has 
driven recent project s such as ATLAS-3D (|Cappellari et all 
1201181 ) and CALIFA jSanchez et al.ll201ll '). These impressive 
projects are limited by their use of monolithic integral field 
units, which means that targeting more than a few hundred 
objects is prohibitively expensive in terms of telescope time. 
This limitation naturally drives us to multi-object integral 
field spectroscopy (i.e. multiplexed IFUs) , the subject of this 
paper. 

2.7 Size and surface brightness 

A key challenge for multi-object integral field spectroscopy is 
to obtain sufficient signal-to-noise at low surface brightness 
levels for all the targets observed simultaneously in a given 
pointing. In this section we present some preliminary inves- 
tigations of the properties of potential SAMI targets. We 
will use the re cent Sersic fits and b ulge-disc decomposition 
carried out bv lSimard et al.l l|201lf ) on the SDSS photome- 
try. 
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Figure 1. The size-surface brightness distribution for an r-band 
sample of SDSS galaxies limited to rgcr < 16.5 (extinction cor- 
rected). Size is defined as the half-light radius, _Re, and the surface 
brightness, /le, is also given at this radius. The marginal distri- 
butions of size and surface brightness are shown below and to 
the right of the main panel. Red dashed lines show the median 
of each parameter, while the dotted red lines show the 10th and 
90th percentiles of the distributions. The blue solid line marks 
the radius of a single fibre core in SAMI and the blue dashed line 
is the radius of the 61-core hexabundle. 



First we consider a simple apparent-magnitude-limited 
sample with rser < 16.5 (extinction corrected, where rscr is 
the galaxy SDSS r-band magnitude derived from a Sersic 
model fit to the photometry). This limit was chosen as the 
galaxy surface density approximately matches the density 
of IFUs in SAMI. The distribution of half-light radius, _Rc, 
and surface brightness at is shown in Figure [1] For such 
a sample the SAMI hexabundle IFUs reach to ITie for all 
but the largest 10 percent of galaxies (blue dashed line), 
while \Ka is sampled by at least 3 IFU elements for all but 
the smallest 10 percent of galaxies. In other words, for the 
central 80 percent of this sample, SAMI can give spatially 
resolved spectroscopy out to at least li?o. The median sur- 
face brightness at ITie is — 22 mag arcsec^ . 

We show how size and surface brightness vary with 
redshift for this sample in Figures [2^ and b. The typical 
sizes in arcseconds of the galaxies stay relatively constant 
with redshift, largely because the r-band limit selects more 
massive (and therefore larger) galaxies at higher redshifts. 
A natural alternative is to select a volume-limited sample, 
which is shown in Figures [2J; and d. In this case we choose 
M-c < —19.5 and z < 0.075, which gives similar numbers of 
targets (i.e. similar surface density) to the r-band cut used 
above. In this case we see, unsurprisingly, that the galaxies 
are smaller at high redshift, but that the median Rc is more 
than twice the radius of a fibre core (and more than three 
fibre cores out to z ~ 0.06). 

An apparent-magnitude-limited sample can be consid- 
ered as a set of volume-limited samples with small redshift 
intervals and absolute magnitude limits that vary with red- 
shift. It is highly likely that the optimal solution for target- 
ing galaxies for a multi-object IFU instrument involves tak- 
ing multiple volume-limited samples, allowing optimal use 
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Figure 2. Size and surface brightness as functions of redshift for 
a magnitude-limited galaxy sample with rsor < 16.5 (panels a 
and b) and for a volume-limited sample with A/r < —19.5 and 
z < 0.075 (panels c and d). The black solid line is the median 
and the black dotted lines are the 10th and 90th percentiles. The 
blue solid line is the radius of a single fibre core in SAMI and the 
blue dashed line is the radius of a 61-core hexabundle. 

of the IFU FoV while broadly populating the distribution of 
galaxy stellar mass. One issue that needs to be considered is 
that in an apparent-magnitude-limited sample, objects in- 
evitably pile up around L* , so suitable sampling may need 
to be implemented to efficiently cover a wide range in stellar 
mass. 

Last considerations in target selection include whether 
to sample all the galaxies in a given volume or to choose field 
locations which uniformly sample the distribution of galaxy 
environment. The challenge with the latter approach is that 
there are a variety of different local density estimators and 
the relation between them is not trivial (see Brough et al., 
in preparation). There is also a challenge to balance the re- 
quirements of maintaining sufficient spatial resolution below 
IRe, as well as reaching > 2i?e which is a requirement for 
reaching the turn-over in the rotation curves of disk galaxies 
(e.g. for TuUy-Fisher analysis). Separate samples within the 
same survey area that can be observed concurrently may be 
the most natural solution to this problem. Further discus- 
sion of the detailed sample selection for a large SAMI galaxy 
survey will be deferred to a future paper. 



3 THE SYDNEY-AAO MULTI-OBJECT 
INTEGRAL FIELD SPECTROGRAPH 
(SAMI) 

With a view to providing the first on-telescope demonstra- 
tion of hexabundle technology, the Australian Astronomi- 
cal Observatory (AAO) and the University of Sydney have 
collaborated on the development of the SAMI instrument 
for the 3.9m Anglo- Australian Telescope (AAT). SAMI uses 
13 X 61-core hexabundles that are mounted on a plug-plate at 
the 1-degree FoV triplet corrector top-end focus of the AAT. 
At f/3.4, with 105 /im core diameter fibres, each hexabundle 
samples a 14.9arcsec diameter field at 1.6arcsec per fibre 
core. At the output end, a total of 13 V-groove slit blocks 
are mo unted at the slit o f the AAT's AAOmega spectro- 
graph l|Sharp et al.l (|2006l ) ; also see Section 13. 7p . Each slit 



block includes 63 fibres (all the fibres from one hexabun- 
dle plus two fibres for sky subtraction). A fusion-spliced 
ribbonised fibre cable of length ~42m joins the two in- 
strument ends together. A near real-time dat a pipeline, 
based on the 2dFdr co de (|Croom. Saunders fc He ald 200j; 
ISharp fc BirchalfcOld ). has been written to reduce the data. 

The following subsections describe the instrument re- 
quirements, fibre cable, the hexabundles, the prime focus 
unit, the SAMI field plates, the AAOmega spectrograph, 
the instrument control and data reduction software. 



3.1 Instrument requirements 

The SAMI instrument was designed to be a technology 
demonstrator and to carry out significant science pro- 
grammes. As a result, the final instrument design is infiu- 
enced by a mix of scientific and technical constraints. A key 
constraint was to develop the system on a rapid time-scale, 
which naturally led to the use of the plug-plate system and 
the already available AAOmega spectrograph. The fiexibil- 
ity of AAOmega, with a range of resolutions and wavelength 
settings also enables a variety of science. 

In order to make a substantial advance over previous 
facilities, the multiplex of the system had to be at least an 
order of magnitude better than previous monolithic IFUs 
(i.e. SAMI required at least ~ 10 IFUs). The median op- 
tical seeing at the AAT is ~ 1.5 arcsec, so the fibre cores 
were approximately matched to this (1.6 arcsec diameter). 
Although undersampling the seeing, this is generally prefer- 
able to having smaller core sizes for a number of reasons: i) 
larger fibre cores provide more independent resolution ele- 
ments, ii) smaller fibre cores will lead to data being read- 
noise limited in the blue, iii) larger fibre cores provide bet- 
ter surface brightness sensitivity, iv) larger fibre cores allow 
higher fill factors, given a fixed minimum fibre cladding (i.e. 
5/im in our case), v) critical sampling of the seeing can still 
be achieved with dithered exposures. 

Given the above, the final key instrument design deci- 
sion was the number of fibres per IFU. The choice of 61 fibres 
per IFU was made largely on the basis of the known capa- 
bility to manufacture such bundles. However, this number of 
fibres also matches other requirements. Larger numbers of 
fibres per bundle would have restricted the multiplex given 
the fixed AAOmega slit length, and the chosen 15 arcsec 
IFU diameter provides a good match to the scale length of 
galaxies in local samples chosen to match the surface den- 
sity of SAMI IFUs within the 1-degree diameter field of 
view (see Section [2.7p . If designed for a single experiment, 
one option would have been to have a variety of bundle sizes 
to match the specific galaxy size distribution of the target 
sample. However, for ease of manufacture, and to maintain 
multiplex and fiexibility, we chose to have all the IFUs the 
same size. 



3.2 Fibre cable 

SAMI is mounted at prime focus on the AAT and feeds 
the AAOmega spectrograph located in the Coude room, re- 
quiring a fibre cable run of ~42m. Within the fibre bun- 
dle a total of 819 fibres are used: 793 for the hexabundles 
(13x61 fibres), and 26 for the sky fibres. The hexabundles 



© 0000 RAS, MNRAS 000, 000-000 



The Sydney-AAO Multi-object IFS 9 



were each supplied with a fibre pigtail length of 1 m encased 
in a reinforced furcation tube. Each of the hexabundle and 
sky fibres were fusion spliced onto a fibre of length ~42 m 
that is terminated in a V-groove block at the spectrograph 
entrance slit (see Section [3.71 below). For this fibre run we 
chose a fibre ribbon cable to minimise the required assembly 
effort. Each 250 micron thick ribbon contains 8 fibres, and 
8 ribbons are used for each of the 13 units (one hexabun- 
dle plus two sky fibres). To minimize losses in the splicing 
process, and because of its ready availability, we matched 
the ribbonised fibre type to that used in the hexabundles 
(ThorLabs AFS105/125Y). 

A 'splice-box' is mounted on the internal wall of the 
telescope top-end barrel that incorporates 13 closed-cell 
polyethylene trays to secure and protect the individual 
splices of each unit. The inputs to this box are the 26 sky 
fibres (each individually sleeved with reinforced PVC fur- 
cation tubing and terminated with SMA connectors) and 
the 13 individual hexabundle tubes. The output of this box 
is the fibre bundle. The bundle is protected during its run 
(which goes through the top end telescope ring, beside the 
primary mirror support, and through the declination drive 
axis to the Coude room) by an inner covering of braided ca- 
ble sleeving and an outer double-split nylon conduit giving 
light-weight yet relatively strong protection. 

Though the use of ribbonised fibre worked well in terms 
of handling and assembly, tests showed that ribbonising 
caused a loss due to focal ratio degradation (FRD), as dis- 
cussed in Section r3.4l The fibre type used provides relatively 
good throughput above 450 nm, though it suffers higher ab- 
sorption losses in the UV than other fibres more commonly 
used in astronomical instruments. The current fibre cable 
described in this paper will be replaced in the first half of 
2012 to provide much improved blue throughput. 

3.3 Hexabundles 

iBland-Hawthorn et al.l (|201ll ) introduced a new imaging fi- 
bre bundle optimized for low-light astronomical applica- 
tions. The hexabundle incorporates several technological in- 
novations in order to achieve a high fill factor. A primary 
goal of the hexabundle design is that the performance of 
individual multimode fibres should be as good as the same 
fibre in isolation. Consequently, we were forced to reject an 
earlier design where the fibre cores were forced into non- 
circular shapes because the added focal ratio degradation 
(numerical apertu re up-conversion) w as found to lead to sig- 
nificant light loss (|Brvant et al.ll201ll ). 

Within each hexabundle, the 61 circular fibres (with nu- 
merical aperture, NA=0.22) are lightly fused together and 
infilled with low-stress glue. The multimode fibres have rel- 
atively small core diameters of 105 /im. In order to achieve 
a high filling factor (75 percent), it was necessary to reduce 
the cladding to only 5 fim. The length of the fused region 
is very short (~30mm) to minimize cross-talk between the 
fibre cores. The fibre bundle is held within a reinforced fiex- 
ible plastic tubing that is both strong and light-weight. The 
optical head is supported by a stress-relieving sleeve and is 
inserted into an SMA connector. This is to allow the fibre 
bundle to be manually attached to the field plate with rela- 
tive ease. The repeatable positioning accuracy made possible 
by the SMA connector is much better than a fibre core size. 




Figure 3. An image of a hexabundle front facet. This shows the 
61 optical fibre cores, which have a total active area of diameter 
980 fj,m. Surrounding the fibre cores is a glass ferrule (black ring), 
which in turn is surrounded by the central steel pin of the SMA 
connection, which extends beyond the edge of the image. 



The use of bare fibre hexabundles directly at the tele- 
scope focal plane is the most novel element of our instrument 
design. There are several advantages to such an arrange- 
ment as an alternative to fibre-lcnslot coupled IFUs that 
have been used elsewhere (e.g. iKelz fc RothI l2006[ l. First, 
the hexabundle is fabricated as a one-step process, whereas 
a fibre-lenslet system requires the fibre bundle to be manu- 
factured and then bonded to the lenslet array (s). Secondly, 
there are no optical elements required in our hexabundle 
design. Lenslet or microlens array IFUs require at least 2 
surfaces and more often up to 8 surfaces before the fibre 
face, to adjust for plate scale and to preserve telecentricity 
of the telescope beam into the fibre. These extra surfaces 
add to system loss. A trade-off must be made against the 
potential for extra losses in the hexabundle (e.g. from FRD - 
see below) and from its lower fill factor. For SAMI, as we are 
feeding the hexabundles with a fast beam the FRD should be 
minimised (if used with the appropriate fibre type) . Finally, 
the hexabundle solution offers the opportunity of reduced 
pitch between IFUs relative to a lenslet solution that must 
work at a higher magnification. 



3.4 Focal ratio degradation and throughput 

Focal ratio degradation (FRD) increases the size of the light 
cone coming out of a fibre compared to that put in. The 
worse the FRD, the more light will be lost from the f/3.15 
acceptance cone of the AAOmega spectrograph. However, 
FRD can be partially controlled by minimising the stresses 
on the fibres when installing them in the system. An addi- 
tional loss comes from the throughput of the fibres, which 
is primarily dependent on the fibre type and length. In this 
section we determine the performance of the fibre cable and 
hexabundles by analysing FRD and throughput. As will be 
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Table 1. FRD and throughput results for the SAMI fibres and 
hexabundles. The bare and ribbonised fibre used is the same AFS 
fibre as in the hexabundles. Results are shown for the central core 
(core 1) of hexabundle number 15 before the hexabundle was 
spliced to the fibre run. Then for core 1 and for 2 other cores (6 
and 18) in the same hexabundle, when the bundle was first spliced 
to the 42 m of ribbonised cable with slit block attached. Lastly, 
for the hexabundle with 42 m of ribbonised cable and slit block, 
but measured after the first SAMI commissioning run, when the 
ribbonised cables had been packed into a braided cable sleeve 
and outer nylon conduit. Cores 1 and 18 are on the edges of a 
ribbonised cable, while core 6 is in the centre of a ribbon. All 
measurements are based on an f/3.4 input to the hexabundle and 
f/3.15 output (accepted by AAOmega). The NA up-conversion is 
the difference between the input and output NA at 90 percent 
encircled energy and has an error in each case of ±0.009. The 
predicted performance is what we expect to achieve when the 
fibre run has been replaced with a higher throughput fibre (e.g. 
Polymicro FBP) and the effects of ribbonising are removed. 



Fibre 


Throughput 


NA up- 


conversion 




Red 


Blue 


Red 


Blue 




% 


% 


90% EE 


90% EE 


Hexabundle 15: 










core 1 alone 


96 ±6 


89 ±5 


0.003 


0.009 


core 1 plus ribbon 


67 ±7 


41 ±6 


0.028 


0.037 


core 1 plus ribbon 










after run 


44 ± 7 


24 ±6 


0.067 


0.073 


core 18 plus ribbon 


64 ±8 


43 ± 7 


0.026 


0.034 


core 18 plus ribbon 










after run 


37 ±8 


22 ± 7 


0.076 


0.082 


core 6 plus ribbon 


65 ±8 


40 ± 7 


0.028 


0.037 


core 6 plus ribbon 










after run 


56 ±8 


32 ± 7 


0.037 


0.044 


Predicted performance 


84 


62 






with replacement fibre 










42m bare fibre 


78 ± 5 


58 ±5 


0.008 


0.018 


42m ribbonised 


73 ±5 


53 ±5 


0.022 


0.031 



seen below, the dominant source of losses is the fibre cable, 
rather than the hexabundles. 

The FRD and throughput were measured for a SAMI 
hexabundle using an LED source that was fed through 
Bessel B and R filters (centred at approximately 435 nm 
and 625 nm) and re-imaged to form an f/3.4 beam. This 
was then input into several hexabundle cores in turn. The 
output fibres were imaged in the far-field using an SBIG 
camera. After fiat-fielding the images, the centre of each 
output spot was fitted. Encircled energy was calculated in 
concentric circles about the centre position using software 
packages within iraf. 

We initially compared the performance of a hexabundle 
alone to the same hexabundle when spliced to the 42 m of 
ribbonised cable with a slit block attached (hexabundle 15, 
'alone' vs 'plus ribbon' in Table [l]). For three different cores 
(numbers 1, 6 and 18; see Figure[3|, the throughput dropped 
to ~ 65 percent in the red and 41 percent in the blue with the 
addition of the ribbonised cable plus slit block. Figure|4]and 
Table [1] show that this drop is identified with a significant 




Numerical Aperture (Input F/3.4 (95%)) 

Figure 4. Encircled energy (EE) vs numerical aperture (NA) 
profiles in B-band for an input of f/3.4. Larger FRD (or NA up- 
conversion) shifts the curves to the right. The vertical line marks 
an output of f/3.15 (into AAOmega). The black curve is the input 
light curve. The cyan line is the FRD of the hexabundle alone, 
without the ribbonised cable. The green, red and blue solid lines 
that are nearly on top of each other, are the curves for cores 1, 
6 and 18 once spliced to the ribbonised 42 m cable and attached 
to the slit block, but before being put into the cable sleeve and 
outer nylon conduit (errors are ±0.006 in NA). The dashed curves 
(with errors of ±0.0065) are for fibres 6, 1 and 18 (same colours as 
above) , after the ribbonised fibre was put into the cable sleeve and 
outer nylon conduit and had been transported to the telescope, 
installed and used for the first commissioning run. 



increase in FRD, with the NA up-conversion (at 90 percent 
encircled energy) for the hexabundle with ribbonised fibre 
measured to be significantly worse than for the hexabundle 
alone. The end finish on the bundle alone was a cleave, but 
the measurement through the ribbonised fibre and slit block 
had the advantage of a polished end finish which should 
improve the FRD, so the FRD introduced by the ribbonised 
fibre may be a little worse than these numbers indicate. 

In order to differentiate between the effect of the 42 m 
of fibre and the ribbonising, we separately compared bare 
and ribbonised fibre throughputs (Table [T] lower panel). 
This was done by measuring the throughput from 10 m of 
bare fibre of the same type as used in the hexabundles and 
ribbonised fibre (AFS105/125Y). The 10 m was then cut, 
spliced together and measured to give the splice loss. Then 
the splice was cut and an additional 42 m was spliced in 
place. This was then cut out and 42 m of ribbonised cable 
was spliced in instead. The throughput could then be com- 
pared for the bare and ribbonised cable after accounting for 
the initial 10 m. Having the initial fibre in place meant that 
the FRD results for the bare and ribbonised fibre were not 
affected by coupling into the fibre or end finishing effects 
as these remained the same for both tests. Variations in 
throughput due to the different splices have been taken into 
account in the errors. Table[l]shows that while the bare fibre 
results will include end effects, the ribbonising of the fibres 
results in significantly worse FRD. The bulk of the loss in 
throughput is due to the length of fibre (78 percent through- 
put in red and 58 percent in blue), however the FRD from 
ribbonising results in less of that light coming out within 
the f/3.15 acceptance cone of the spectrograph. 
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Figure 5. Schematic of AAT prime focus triplet corrector. 



Once the hexabundles were spliced to the ribbonised ca- 
ble and slit blocks ('plus ribbon' in Table[l|, the throughput 
was lower than that measured for the 42 m of ribbonised ca- 
ble alone. This is due to an increase in the FRD above that of 
the ribbonising. The small additional FRD is therefore due 
to the hexabundle and slit block. Mechanical stress on the 
fibres when moun ted in the slit block can increase the FRD 
(see for example lOliveira. de Oliveira &: dos Santos! l2005h . 
and hence decrease the throughput within f/3.15. 

During the commissioning run it was noted from 
AAOmega images that the throughput of the fibres on the 
outer edge of each ribbonised cable was lower than that of 
the fibres in the centre of a ribbonised cable. Therefore two 
fibres on the edge of a ribbon and one in the centre of a 
ribbon were tested after SAMI had been taken off the tele- 
scope. The testing method was identical to that described 
above, and the three fibres tested were the among the same 
fibres from the bundle that was tested before the SAMI com- 
missioning run. In Figure [4] and Table [T] it is notable that 
while the fibres had similar performance before commission- 
ing (see 'core plus ribbon' in the Table), after commissioning 
(see 'core plus ribbon after run') the FRD was substantially 
worse for the two fibres on the edge of a ribbon (fibres 1 and 
18). Meanwhile, the fibre in the centre of the ribbon (fibre 
6) showed only a small increase in NA up-conversion com- 
pared to the previous results. The throughput within f/3.15 
for the cores on the edge of a ribbon decreased from 64-67 
percent and 41-43 percent (red and blue respectively) before 
the SAMI installation to 37-44 percent and 22-24 percent af- 
ter the SAMI commissioning run. However, the hexabundle 
core in the centre of a ribbon (core 6) had a comparatively 
higher throughput of 56 percent and 32 percent. 

We believe that the increased FRD for the edge fibres 
is due to stresses in the ribbonised cable when it was fed 
into the braided cable sleeving and outer nylon conduit that 
protects the fibre run. As the fibre bundle was moved during 
transport, installation and commissioning, the ribbonised 
bands could be bent in all directions, in which case the edge 
fibres come under more stress than the centre. 

The throughput of SAMI will soon be significantly im- 



proved by the replacement of the fibre run. The loss due 
to the fibre length will be reduced, particularly in the blue, 
by using a different type of fibre. The fibres being consid- 
ered will result in up to 14 percent higher throughput at 
the blue end of the spectrum. In addition, if the effects of 
ribbonising and FRD stresses from the ribbonised cable are 
removed, the throughput would further increase by up to 
~27 percent, perhaps doubling the current blue through- 
put. An alternative protection will be required for the 42 m 
fibre run, and that may introduce some FRD losses, but it is 
being designed to have less of an effect than the ribbonising. 

3.5 Prime focus unit 

Within the AAT prime-focus top-end, the triplet corrector 
and the Prime Focus Camera were refurbished for the SAMI 
instrument. Originally built for the commissioning of the 
AAT in the early 1970s, they provide a 1-degree FoV with a 
plate scale of 15.2 arcsec/mm. With 61x105 /im cores, this 
provides a 15 arcsec diameter FoV for each hexabundle and 
a sampling of 1.6 arcsec per fibre core (see Figure [5)|. 

With a total of 39 fibre positions in each field (26 sky 
positions and 13 object positions), we chose to use a plug- 
plate assembly rather than a robotic positioning system, as 
the operational overhead and down-time between fields for 
reconfiguration are both relatively low, particularly as we 
are targeting long (~2-hour) integrations and have a con- 
nectorised fibre system. 

The SAMI plug plates are pre-drilled 3 mm-thick brass 
discs with through- holes at each object/sky location. Each 
hexabundle and sky fibre is terminated in a SMA screw- 
thread fibre connector. A mating connector is installed in 
each plug-plate at each position. Two galaxy fields (i.e. 26 
objects) are pre-drilled on the plate along with a set of 26 
blank-sky locations common to both galaxy fields. For the 
proposed integration time of 2 hours per field, this means 
that 2 plates (and 1 plate exchange) are required each ob- 
serving night. The down-time between fields recorded during 
the initial commissioning run was less than 30 minutes. Fur- 
ther investigation found that it was possible to include 4 
fields per plate, so that no plate exchange is required during 
the night. 

The plug-plates are installed within an assembly that 
is kinematically mounted to the Prime Focus Camera (see 
Fig. E}. Due to the relatively large FoV of each hexabundle, 
there is not a strict requirement on the positioning accuracy 
of each hexabundle central fibre. However, we aim to reach a 
total positional accuracy of half the core diameter of a single 
fibre (i.e. 0.8 arcsec). The accuracy achieved is determined 
by several factors that include: the hexabundle concentricity; 
the connector concentricity; the plug-plate machining accu- 
racy; and the plug-plate thermal expansion. These factors 
are all controllable to within much less than a fibre core. 
Additional position errors are introduced via i) a rotation 
offset, which is corrected via a fine-thread micrometer ro- 
tation adjustment between the plug-plate assembly and the 
Prime Focus Camera, ii) an asymmetric radial error arising 
from an x-y offset of the plug plate from the telescope op- 
tical axis, for which the tolerances are quite large and iii) a 
symmetric radial error arising from an incorrect plate scale, 
which can be corrected after measurement by producing a 
new distortion map. 
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Figure 6. The SAMI plug-plate assembly unit mounted onto the 
Prime Focus Camera. The white 'splice box' (top) connects the 
blue hexabundles and orange sky fibres from the brass plug-plate 
to the fibre bundle. The cross-bar above the field plate provides 
a location to mount the guide camera, which images the hole in 
the middle of the plate. 

For image reconstruction it is important that the rota- 
tional alignment of each hexabundle is known. Because the 
SMA connector has no rotation adjustment capability and 
is not rotationally keyed, the orientation of the hexabundles 
was determined by eye as they were inserted into the plug 
plate. Laboratory and on-sky tests demonstrated that an ac- 
curacy of less than half a core in the outer hexabundle ring 
(<10 degrees) was possible with this technique. The SMA 
connectors will shortly be replaced by keyed FC connectors 
to eliminate this limitation. 

For acquisition and guiding, we use a CCD camera 
mounted on a gantry above the plug-plate that views (via an 
optical relay) the central region of the field through a hole 
in the plate. This camera (an 800x600 pixel Watec #120N) 
provides a video output signal that can be integrated from 
frame rate up to ~10 seconds and is compatible with the 
existing AAT control system guiding software. It reaches 
~14th magnitude over a sky FoV of ISOarcsec diameter 
(providing a sufficiently high sky coverage factor for most 
galactic latitudes) and has a sampling of 0.3 arcsec/pixel. 



3.6 Field plate manufacture 

Using the manual plug-plate method for positioning the hex- 
abundles in SAMI presents a number of challenges for the 
field allocation and plate manufacturing. Here we describe 
the SAMI field configuration methodology. 

The plates are manufactured using brass of thickness 
3 mm and approximate diameter 240 mm (corresponding 



to the 1 degree FoV). The fibre connectors have an as- 
signed footprint of 15 mm to allow access for installation 
and removal. The instrument uses two distinct physical plate 
types: science and calibration. Each plate is configured with 
multiple stacked fields for efficiency gains. There is a central 
10 mm hole in each plate for the guide camera to image a 
sky region of diameter ~150 arcsec. 

The science plates consist of two stacked galaxy fields 
each with 13 galaxy targets for the IFUs and about 5 field 
alignment stars, also targeted with IFUs. The field align- 
ment stars did not share the same field centre as the galaxy 
targets. During commissioning, when the plate was first ob- 
served, IFUs were positioned at the location of the align- 
ment stars in order to check for consistent rotation and plate 
scale between each plate. Once confirmed, the IFUs were re- 
allocated to galaxy targets and the telescope moved to the 
science field centre. Each plate also contained 26 sky fibre 
positions shared between the two galaxy fields. The science 
plate is configured for consecutive fields requiring only the 
re-positioning of the 13 IFUs onto new galaxy targets. Each 
science plate then has a total of ~ 62 drilled holes for the in- 
stallation of the SMA fibre connectors (with a minimum sep- 
aration of 15 mm). The astrometric calibration plates were 
constructed in a similar way, but in that case each plate 
contained four overlapping fields and four visual alignment 
fields. 

The process for generating a single field was as follows: 

(i) For each RA region, convert the RA and Dec. coordinates 
to angular distance coordinates and calculate the pairwise 
distance matrix of all targets. 

(ii) Iterate over each target using the distance matrix to ex- 
tract targets within a 0.5 degree radius and to remove targets 
with separations <15mm (~228 arcsec). 

(iii) Count the number of targets and identify this as a can- 
didate field if the total is equal to or greater than A'', where 
A*' is the number of sources required; A'^=14 for astrometric 
calibration fields (13 IFU targets plus central guide star), 
N=3 for visual alignment fields (two visual targets and a 
central guide star) and A''=13 for science fields. 

(iv) For a candidate field, save the data and assign the field 
to that RA region; eliminate the field's allocated objects 
from the list of potential targets (in order to produce fields 
with unique sets of targets). 

Once individual fields were defined, multiple fields were 
stacked on a single plate (e.g. for science fields the plates 
contain 2 x 13 galaxies plus 2x5 alignment fields; for as- 
trometric calibration fields the plates contain 4 x 13 calibra- 
tion targets plus 4x2 visual alignment fields). Fields were 
grouped into RA regions (typically of size 20 degrees in RA) 
based on observation times. The RA regions were ordered 
based on the best observing time for each region (i.e. the 
time of lowest airmass). Then fields were tested in turn to 
see that the allocated targets did not overlap, so that they 
could be stacked onto an individual plate. The priorities of 
the targets in each science field are summed to give a figure 
of merit for all possible plate permutations. The permuta- 
tions can be filtered for field duplicates and sorted by merit. 

Once the plate configuration was defined, we deter- 
mined the 26 sky fibre positions that are shared for each 
field. This was based on a grid of regularly spaced angular 
positions relative to the plate centre (e.g. 25 x 25 grid points 
over the field), which were then converted to RA and Dec. 
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SAMI Plate - comm1_galaxy_p4.dat 

View looking at Plate from PF access position. 




Figure 7. Example SAMI science plate schematic showing 2x13 
galaxy fields (yellow/orange), 26 common sky fibres (blue), 2x5 
field alignment stars (pink) and the central guide window (green) . 



given each field centre, and to relative plate coordinates in 
microns (the same for each field). Sky fibre grid locations 
were eliminated if they overlapped with the science targets 
in pla te coordinates. SuperCOSMOS images (|Hamblv et al.l 
were examined at the RA and Dec. of each remain- 
ing sky grid position for each field and if no source was 
found within a 24arcsec window then the sky position was 
accepted. From these remaining candidates, 26 sky fibre lo- 
cations were then selected for targeting. 

For each stacked field, we then determined suitable 
guide stars {V<14) in the central 150arcsec FoV using the 
Aladin sky atleis | http: / / aladin.u-strasbg.fr/| to search the 
USNO-Bl catalogue and examine DSS images. 

After the above process, we take the entire set of tar- 
get positions in RA and Dec. and convert them to accurate 
plate X-Y coordinates using the optical distortion model of 
the prime focus corrector and the differential atmospheric 
refraction at the expected time of observation. These X-Y 
positions were then converted to mechanical drawings (with 
temperature compensation to allow for the difference be- 
tween the average day-time drilling temperature and the 
average night-time observing temperature), and these were 
the input to a CNC machine for drilling the plates to the 
required specification. A schematic of a galaxy field plate is 
shown in Fig. m 



3.7 AAOmega spectrograph 

AAOmega (jSharp et al.l |2006| ) is a fibre spectrograph de- 
signed for use with the 2-degree Field (2dF) robotic fibre 
positioner on the AAT. It has a double-beam Schmidt design 
that allows for optimised performance in the red and blue 
simultaneously. Volume-phase holographic (VPIS) gratings 
are employed to give high throughput. There are a range of 
gratings available, giving resolutions from ~1700 to '^13000 



for the SAMI fibres (105 fim diameter). There is a choice of 
two dichroics that split the light and direct it to the blue 
and red cameras; one dichroic has a cut-off wavelength of 
570 nm and a second at 670 nm. The detectors are 2kx4k 
E2V devices with the short axis in the wavelength direction 
and the long axis in the spatial direction. 

For operation with the 2dF top-end fibre feed, the 
AAOmega slit is populated with 392 fibres that each have 
140 /im cores, projecting to ~3. 4 pixels on the detector, with 
a core-to-core pitch of 10 pixels. AAOmega is also fed by the 
SPIRAL IFU used at the AAT Cassegrain focus. SPIRAL 
has a 16x32 rectangular lenslet array, with a sampling of 
0.7arcsec. The 512 SPIRAL fibres have 85 /xm core diame- 
ters, which project to ~2. 4 pixels on the detector. The fibre 
cores are separated on the slit by a pitch of 140 /im. 

For SAMI we have matched the fibre core-to-pitch ra- 
tio of SPIRAL, which is 85/140 = 0.61. This is more 
closely packed than used for 2dF, but has proven sufficient 
with SPIRAL for the minimisation of crosstalk between fi- 
bres using an opt i mal ex traction methodology described by 
ISharp fc Birchalll (|2010[ ). The SAMI fibre pitch on the slit 
is 170 /im, giving a core-to-pitch ratio of 105/170 = 0.62. 
This pitch then defines the maximum number of hexabun- 
dles (13) that will fit along the concave AAOmega slit while 
allowing 26 individual sky fibres to also be positioned on the 
slit. 

For convenience of manufacture and assembly, each hex- 
abundle and sky fibre pair is fed to an individual silica V- 
groove slitlet. The mapping is such that the central core is in 
the middle of the slit. The outward spiral numbering of the 
hexabundle cores corresponds to odd and even outwardly 
alternating fibres in the slitlet, with the sky fibres being the 
outermost fibres on the slitlet. This arrangement maximises 
the number of adjacent hexabundle cores that are also ad- 
jacent on the detector and minimises cross-contamination 
between sky and adjacent object spectra. 

The AAOmega slit is mounted on a rotating mechanism 
that allows four different slits to be positioned in the beam 
of the spectrograph. This is used during operation of 2dF, 
when one of the two fields is being observed while the other is 
being re-configured by the robot. The SAMI slit is located 
at one of the previously spare positions on the AAOmega 
slit rotator; in the course of normal SAMI operations the 
slit does not move. 



3.8 Control software 

SAMI uses an adaption of the 2dF/ AAOmega control soft- 
ware (Smith et al. 2004). This software has significant flex- 
ibility and had been successfully evolved through various 
changes since it was originally commissioned with 2dF in 
1996. As part of the implementation of AAOmega in 2006, 
the software was modifled to support the SPIRAL IFU feed 
to AAOmega. This operational mode did not use the flbre 
positioner and ADC of the 2dF top-end, and as a result 
was a good starting point for SAMI. A number of minor 
modifications were required to the software to support this 
additional mode of operation. 

The 2dF/ AAOmega software inserts a large binary ta- 
ble into the CCD data files to fully describe the location 
and status of each of the fibres. This table differs signif- 
icantly from versions implemented for SPIRAL and 2dF. 
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A component of the software (the Fibre To Fits program) 
is used to build this table, and a new version of this pro- 
gram was required for SAMI. It reads a file describing the 
allocation of objects to IFUs and sky fibres. The software 
combines this with the telescope pointing information and 
information from files containing the measured positions of 
each fibre in each IFU to determine the actual location of 
each fibre on the sky. A s a result, a EURO-3D compliant 
l|Kissler-Patig et al.ll2004l ) binary table can be generated and 
World Coordinate System (WCS) information given for each 
fibre. Extra information in the table is provided to ensure 
every fibre is fully traceable, and that its position and the 
origins of its position are well defined. 

It is important that fields be observed sufficiently close 
to the time planned, and for which the plug-plate was drilled. 
Otherwise changes in airmass can cause some objects to be 
poorly acquired, largely due to differential atmospheric re- 
fraction across the FoV. To assist in the decision whether 
to observe a given plate, the control software displays the 
error between the drilled hole and the current object sky 
position for each of the 13 IFU probes at the beginning of 
each exposure. 



3.9 Data reduction software 

Data reduction for SAM I is performed using the 2dFdr 
data reduction pipelin e (|Croom. Saunders fc Healdl 12004 
ISharp fc BirchaOlbOlOl ) originally written for the 2dF instru- 
ment and then modified for use with the AAOmega spectro- 
graph. This provides fully automated reduction of flat fields, 
arcs and object frames, including spectral extraction, wave- 
length calibration and sky subtraction. The code is run from 
a fully configurable GUI which allows user control of the al- 
gorithms used. The main modifications required for SAMI 
were a revised algorithm to accurately map the fibre loca- 
tions across the detector (the 'tramline map') and routines 
to read new elements of the FITS binary table in the data 
frames. 2dFdr is sufficiently fast that data can be reduced in 
real time, a feature that is of particular value during com- 
missioning of an instrument. 

The 2dFdr pipeline generates extracted, flat-fielded, 
wavelength-calibrated, throughput-normalized and sky- 
subtracted spectra. The data product is a 2D image con- 
taining each of the fully reduced ID spectra, together with 
a variance array and binary table (see Section FS.SI) . To allow 
fast reconstruction of 3D data-cubes and first-pass science 
analysis, we have developed a suite of python-based routines 
which allow us to: (i) construct and write full data cubes for 
individual IFUs (ii) view all 13 reconstructed IFU images 
collapsed over a user-specified wavelength range; (iii) calcu- 
late the centroid of a source within the IFU field; (iv) cal- 
culate the offset in arcseconds between the source and the 
centre of the IFU; (v) extract summed spectra for an entire 
IFU or view individual spaxel spectra; and (vi) fit emission 
lines and construct kinematic maps on the fiy. 

As an aid to the development of the data reduc- 
tion pipeline, an instrument data simulator was used. This 
was based on the simulator used for the new high resolu- 
tion HERMES spectr ograph being developed for the AAT 
ijGoodwin et al.ll201Gl ). but modified to simulate SAMI de- 
tector images. This was of particular value in developing the 



algorithm required to extract the fibre spectra from the 2D 
image. 

3.10 Hardware upgrades 

The current system, as described above, performs well (see 
results below), but substantial improvements can be made 
in performance and usability. To this end a new fibre cable 
will replace the current one in the first half of 2012. This 
will substantially improve the system throughput, particu- 
larly at wavelengths < 4500A. The new fibres will be housed 
in a low stress cable to further improve throughput by sub- 
stantially reducing FRD. A third planned modification will 
be to use a different form of connector for the hexabundles. 
While providing accurate alignment in the x and y, the cur- 
rent SMA connectors do not allow easy rotational alignment. 
The SMA connectors will be replaced with a keyed connec- 
tor such as the FC type. 



4 COMMISSIONING OBSERVATIONS 

4.1 Stellar target selection 

The astrometric calibration plates consist of four stacked 
stellar fields (Vt < 10 mag), each having 13 IFUs for dis- 
tortion mapping and alignment refinement as well as four 
stacked bright stellar fields ( Vt < 6.5 mag) each with at 
least two stars for initial visual alignment. The calibration 
plate field uses a centred guide star located at the field cen- 
tre. 

The calibration plates are derived from the stellar data 
contained withi n the Tycho-I Reference Catalogue (TRC; 
iHog et al] 1 19981 ). We trimmed the catalogue on position 
(Dec. < 20 degrees) and magnitude (Vt < 6.5 for visual 
alignment fields, Vt < 10 for calibration fields) for a num- 
ber of observable RA regions. 

4.2 Galaxy target selection 

The targets for the commissioning were selected to test the 
full capability of SAMI, subject to the following criteria: 

(i) The density of targets on the sky exceeded the 13 hex- 
abundle units available in the 1-degree field of SAMI. 

(ii) The targets were of sufficient angular size to fill the 15 
arcsec aperture of each fibre bundle. 

(iii) Galaxy surface brightnesses (in either emission lines or 
continuum) were likely to yield S/N>3 in the outermost fi- 
bres of each bundle in a ~2-hour integration. 

We decided to select commissioning targets from a wide 
cross-section of galaxy spectral types to better constrain the 
performance of the instrument across a range of future sci- 
entific applications. 

The target input catalog ue was drawn from the 6 - degre e 
Field Galaxy Survey (6dFGS: [jones et al.|[2003 . l2005l . 120091 ') . 
which itself was selected from the Two-Micron All-Sky Sur- 
vey E xtended Source Catalog (2MASS XSC; I Jarrett et al] 
'2000!). The 6dFGS catalogue contains 125071 redshifted 
southern galaxies (|6| > 10°) selected to K < 12.65, H < 
12.95, J < 13.75, rp < 15.60, and 6,7 < 16.75. The median 
redshift of the sample is z = 0.05. Apparent galaxy sizes are 
taken from the 2MASS XSC. 
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Table 2. Priorities used in allocating targets for the SAMI commissioning observations (9 is the highest priority and 1 is the lowest). 



Priority 


Input sample 


Spectral criterion 


Radial size range 


Targets 


9 


6dFGS main 


strong emission lines 


r > 10" 


100 


8 


6dFGS main 


strong emission lines 


5 < r < 10" 


748 


7 


6dFGS-v 


high-S/N early-type 


r > 10" 


28 


6 


6dFGS-v 


high-S/N early-type 


5 < r < 10" 


152 


5 


6dFGS main 


strong emission lines 


4.46 <r <5" 


107 


4 


6dFGS-v 


high-S/N early-type 


4.46 < r < 5" 


202 


3 


6dFGS-v 


high-S/N early-type 


r < 4.46" 


216 


2 


6dFGS main 


strong emission lines 


r < 5" 


298 


1 


6dFGS main 


all others remaining 


all others remaining 


2235 








Total 


4086 



The mean density of 6dFGS galaxies is ~7deg ^, sub- 
stantially lower than the ~19deg"^ IFU density of SAMI. 
However, the low-redshift nature of the sample means that 
6dFGS target densities vary considerably over the sky, and 
in the densest regions, 6dFGS target densities exceed the 
IFU density of SAMI. By careful pre-selection of dense re- 
gions we ensured that all 13 hexabundles were filled for every 
field. Nine dense regions (each of diameter 8 deg) were se- 
lected over a range of hour angle and declination to provide 
4086 potential targets for the final allocation of SAMI fields 
and IFUs. 

Targets were ranked on a scale from 1 (lowest priority) 
to 9 (highest priority) that was used to weight the targets 
for hexabundle allocation and field placement. Targets were 
given a greater relative weighting if they had: (i) large an- 
gular size; (ii) prominent spectral emission lines; or (iii) a 
spectrum typical of an early-type galaxy, with high signal- 
to-noise. 

The angular sizes used in criterion (i) were 2MASS J- 
band half-light radii. Galaxies were divided into those with 
apparent radial sizes r > 10" (logr > 1.0), those with 
5 < r < 10" (0.7 < logr < 1.0), those with 4.46 < r < 5" 
(0.65 < logr < 0.7), and the smallest ones that remained, 
with r < 4.46" (logr < 0.65). These divisions were chosen 
to differentiate between those galaxies most closely matched 
to the hexabundle aperture size (with galaxy diameters 10- 
20"), and those outside this range. For the purpose of com- 
missioning, the largest galaxies were also given high priority 
to extend the range of test subjects. 

Criterion (ii) was based on matches to the emission- 
line galaxy spectral templates used duri ng the 6dFGS cross - 
correlation redshifting procedure (see Ijones et al. 1 l2004h . 
Note that it is not a complete sample, as spiral galaxies 
which are dominated by their bulge on the scale of the 
6dFGS fibre aperture (6.7" diameter) may not exhibit sub- 
stantial emission lines in their 6dF spectrum. 

Criterion (iii) was populated by the 6dFGS velocity c at- 
alogue (6dFGS-v: ICampbelll [2OO9I : ISpringob et all 120111 ). a 
high signal-to-noise (> 10 per pixel in the 6dFGS spectrum) 
subsample of 11288 early-type galaxies. 

Selecting the target sample in this way ensured that 
the performance of SAMI would be gauged across a broad 
range of galaxy types, while also maximising size and sur- 
face brightness considerations, where possible. Table [2] sum- 
marises the selection criteria applied to each of the priority 



assignments, and the total number of available targets in 
each case. 

4.3 Observations 

The first commissioning observations were carried out on the 
nights of 1-4 July 2011 at the AAT. The primary aims were 
to test the alignment and astrometry of the plug-plate and 
hexabundle units, estimate system throughput, and exam- 
ine the reliability and robustness of the instrument. The sec- 
ondary goal was to obtain the first galaxy IFU observations 
with SAMI to test data quality and determine the accuracy 
with which physical parameters could be extracted. 

The AAOmega spectrograph setup used the 580V grat- 
ing in the blue arm at a central wavelength of 480 nm, cover- 
ing the wavelength range from 370 nm to 570 nm (the latter 
set by dichroic that splits the light between the red and blue 
arms of AAOmega). With the 105 /im fibre cores of SAMI, 
the 580V grating provides a spectral resolution of J? = 1730 
(or 173 kms"^ FWHM). The dispersion is 0.103 nmpixeP^ 
This wavelength range is ideal for measuring a wide range 
of spectral features in the blue parts of galaxy spectra at 
low redshifts (z ~ 0.05), such as the D4000 break, [O III] 
emission, various hydrogen Balmer lines, and Mgb. In the 
red arm, the key spectral features are the Ha, [N II] and 
[S II] emission lines. As these are all located in a relatively 
narrow wavelength range (and the target redshift range is 
small), it is possible to observe all of these features using 
the higher resolution lOOOR grating. This was set to have 
a central wavelength of 680 nm, providing a spectral range 
from 625 nm to 735 nm. The spectral resolution in the red 
arm is R — 4500 (or 67kms~^ FWHM), with a dispersion 
of 0.057 nmpixeP^. 

The first observations used astrometric calibration 
plates to measure the accuracy of IFU placement, including 
checks of rotation, plate-scale and distortion. An initial cor- 
rection for rotation was applied by one of us (SR) observing 
a small number of bright (Vr < 6.5) stars by eye from within 
the AAT prime focus unit, using calibration holes drilled in 
the plate at the position of the stars. The stars were aligned 
to the holes using the rotational adjustment micrometer on 
the plate holder assembly (Section 13. 5|) . Once this first cor- 
rection was made, several astrometric calibration fields were 
observed in order to make precise measurements of the re- 
quired parameters. 

After the astrometric observations where carried out. 
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Figure 8. The relative fibre-to-fibre throughput of the hexabun- 
dle (crosses) and sky (red circles) fibres as a function of fibre num- 
ber. Damaged hexabundles are not displayed (including 3 where 
the damage is so great the data is not useful, and a 4th where al- 
though the throughput is low, there is still reasonable S/N in the 
fibres). Inset is a small part of the throughput distribution from 
fibres 70 to 120. Here the periodic decline in throughput at the 
edges of each fibre ribbon can be seen. The dashed vertical lines 
break up the fibres into separate fibre ribbons, each containing 8 
fibres. 

we targeted a number of standard stars to measure system 
throughput, and then a galaxy field was observed as part 
of science verification. During periods of poor weather other 
tests were carried out, including examination of guide cam- 
era flexure and variations in fibre throughput as a function of 
telescope position. The detailed results of all commissioning 
observations are presented in Section [5] below. 



5 RESULTS 

5.1 Data reduction 

All data presented below was reduced using 2dFdr (see Sec- 
tion l3.9|l . Fibre flat fleld frames illuminated by a quartz lamp 
are reduced first to define the location of the spectra ('tram- 
line map'), and to construct a fibre-flat fleld, which all other 
frames are normalized by. This provides correction of the 
relative colour response of the fibres, but not total through- 
put normalization. Extraction of the spectra was carried out 
using Gaussian fits to the fibre profiles, based on a mean 
profile shape. Wavelength calibration was via an arc frames 
using a copper-argon lamp. Throughput calibration of the 
fibres was carried out either by using twilight sky frames or 
the strength of night sky emission lines in the object spec- 
tra. Once the fibres were flat-fielded, wavelength-calibrated 
and corrected for relative throughput, sky subtraction was 
performed. A median sky spectrum was constructed from 
the 26 sky flbres. This was then subtracted from the object 
spectra. 

The relative throughput of each flbre is shown in Fig|8l 
These are normalized to the median throughput, so are dis- 
tributed about a value of 1. The full range of throughputs is 
from ~ 0.8 to 1.2. A marked periodic structure is seen every 
8 flbres (see inset in Fig[8ll which is due to poorer throughput 
at the edges of each fibre ribbon (see Section [3.4p and the 



relative difference is consistent with the expected through- 
put loss from FRD measured in the lab based experiments. 
The measured rms variation in relative fibre throughput is 
0.13, however the uniform throughput away from the ribbon 
edges suggests that a rebuilt fibre cable not using ribbons 
would provide much more uniform throughput. 

Good throughput calibration and sky subtraction are 
crucial in obtaining accurate spectroscopy at faint surface 
brightnesses. Here we make a preliminary assessment of sky 
subtraction precision, noting that various upgrades to the 
data reduction code are still to be implemented. Here we fo- 
cus on the residual continuum in fibres which were located 
on blank sky. Minimizing the residual continuum is criti- 
cal in enabling accurate measurements of stellar population 
parameters (e.g. age and metallicity) from absorption line 
indices. The subtraction of night sky emission lines, while 
also important, can be tackled in a va riety of ways, includ- 
ing principle component analysis fiSharp fc Parkinson 201(3) 
that can effectively remove emission to the Poisson limit 
of the data. This approach will be implemented for SAMI 
data. As the galaxy commissioning targets completely filled 
the field of view of each hexabundle (see Figs. llll and ll6|l . we 
used only the sky fibres for this test. It should be noted that 
examination of the spectral point-spread function (PSF) of 
sky and hexabundle fibres found no noticeable difference be- 
tween them, suggesting that a test of only the sky fibres pro- 
vides a fair assessment of sky subtraction accuracy. Future 
on-sky tests will confirm this with blank sky observations 
using the hexabundles. We median filter the sky spectra be- 
fore and after sky subtraction and then take the ratio of 
the sum of these median filtered spectra as a measure of 
the fractional residual sky continuum present. The rms con- 
tinuum residual about zero is 3.5 percent, suggesting that 
in the current data this is the level at which we are accu- 
rately subtracting sky continuum. There are a small number 
(typically ~ 2) sky fibres which perform significantly worse 
that this (residuals of ~ 10 percent or more). Examination 
of these fibres showed that i) they do not demonstrate sky 
emission line residuals at the same level, and ii) they tend 
to be located next to the most badly damaged hexabundle 
which projects essentially no light onto the spectrograph 
CCD. This points to residual scattered light in the spec- 
trograph as the main cause. In particular, the fibre profiles 
are made up of a Gaussian core and low level, but broad, 
scattering wings. In regions of the CCD which are fully pop- 
ulated with fibres the scattering wings coadd to form an 
approximately constant pedestal above the bias level of the 
detector (when data are sky limited). However, in detector 
regions without illuminated fibres the background level falls 
below the pedestal caused by the scattering wings. As a re- 
sult, fibres on the edge of a blank region have a modified 
PSF and can be poorly extracted from the data frame. The 
solution to this is to fit a m ore complex PSF to th e data 
in the extraction process (e.g. lSharp fc Birchall|[201ol ). This 
will be implemented within the 2dFdr package. 

5.2 Astrometry 

The AAT Prime Focus triplet corrector distortion is dom- 
inated by a pincushion effect. For SAMI, we implemented 
a model of th is using the SLA library routine SLA_PCD 
(| Wallace! [T99i ) and the recommended value of 178.585 for 
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Figure 9. The overall measured throughput for the SAMI sys- 
tem in the blue arm of the AAOmega spectrograph using the 
580V grating. The red star shows the measured throughput av- 
eraged over the band indicated by the green horizontal line. The 
purple star shows the lab test estimate, with error bars. The hor- 
izontal cyan line indicates the wavelength range over which the 
laboratory measurements were made. The low throughput in the 
blue was expected from our choice of fibre for the demonstrator 
instrument. 



the pincushion coefficient, c. The radial distance in the pres- 
ence of distortion is given by 



■ r(l + cr 



(1) 



where r is the radial distance from the tangent point. Dis- 
tances are in units of the projection radius. An analysis of 
a Zeemax model of the optical system agreed with the pre- 
viously used model values, to a maximum error of 0.25 arc- 
seconds at the field edge. Careful attention was paid to the 
many orientation and sign issues throughout the software 
and plate manufacturing process, leading to the first plate 
being successfully acquired with only a minor rotation cor- 
rection required. 

The AAO implemented a package for calibrating the as- 
trometric models in multi-o bject spectrographs fo r the ESO 
VLT FLAMES instrument (iPasguini et ahlbOOj) and later 
for the Subaru FMOS instrument (|Kimura et al.lboiol '). This 
package, known as FPCAL, can be configured for new in- 
struments by the use of plug-in software components. These 
components provide the program with an implementation of 
the optical model and the ability to read details of objects 
allocated to fibres and the telescope astrometric model pa- 
rameters from the FITS files generated by the instrument. 
FPCAL then provides tools to allow fitting of model param- 
eters and analysis of the results. The user of FPCAL can se- 
lect either Singular Value Decomposition (SVD) or Powell's 
method of minimization (e.g. [Pre ss, Fla nnerv fc Teukolskvl 
1 19861 ) as the fitting algorithm. In practice, these give the 
same results within noise margins, so the SVD technique 
is normally chosen as it is faster and provides more infor- 
mation. The appropriate plug-ins to enable us to calibrate 
SAMI with FPCAL were implemented. FPCAL requires as an 
additional input the offset necessary to centre each IFU, and 
software was implemented in python to extract this from the 
reduced data files. Centroids were calculated using a Gaus- 
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Figure 10. The overall measured throughput for the SAMI sys- 
tem in the red arm of the AAOmega spectrograph using the lOOOR 
grating. The green star shows lab test estimate for the throughput 
over the ij-band (5750-7000A) indicated by the blue horizontal 
line. The comparison between directly observed and lab-predicted 
throughputs is approximate as the on-sky observations did not 
cover the whole of the photometric ij-band. 



sian fit to determine the position of calibration stars with 
respect to the IFU centre. 

Observations of three different astrometric calibration 
fields were made, including fields drilled in different physical 
plates. An analysis that includes acquisition error, scale, ro- 
tation and distortion was performed for each individual field. 
Sets that combine multiple observations where the stars were 
offset by ~2 fibre cores with respect to the IFU centres were 
used to test the rotation of individual IFUs. 

In the initial analysis, results were entirely dominated 
by minor acquisition and scale issues; the following refers 
to one typical observation. An acquisition error of approxi- 
mately 1.25" was removed to simplify further analysis. A fit 
to the scale found a focal length about 20 mm longer than 
originally presumed; this scale error caused positional errors 
of up to 1.2" at the edge of the plate. Removing the scale er- 
ror resulted in 0.8" RMS residuals. Fitting the pin-cushion 
distortion parameter and field rotation did not provide a 
significantly better fit (0.78" RMS residuals). Whilst the 
change in focal length of 20mm is larger then expected, this 
is believed to be due to the power of the triplet corrector, 
and the actual physical change involved is likely small. 



5.3 Throughput 

During the SAMI commissioning run two spectrophotomet- 
ric standards were observed over two nights in several of the 
individual IFUs. On 2 July 2011 the star LTT6248 was ob- 
served in IFU hexabundles 10 and 13 (H#010 and H#013). 
The first of these observations was used to calculate the 
overall throughput of the SAMI system including the at- 
mosphere, telescope primary, prime focus corrector, SAMI 
fibre feed, and AAOmega spectrograph. First an integrated 
spectrum was extracted from the IFU data by summing the 
spectra in all spatial elements of the IFU. Since the data 
are under-sampled and the IFU grid is not contiguous, some 
light is lost between spaxels and a summed spectrum can 
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Table 3. Galaxies observed with SAMI during the July commissioning run. Uncertainties on magnitudes and colours are typically ~ 0.1 
mags. 



Hexabundle 


R.A. 


Dec. 


Priority 


6dFGS ID 


Redshift 


bj 


bj - rp 


bj-K 


Radius^ 


Throughput 




(J2000) 


{J2000) 














(arcsec) 


n 9 


H#016 


19 55 18.8 


-54 59 35 


1 


gl955188-545935 


0.019 


15.43 


1.29 


4.68 


11.6 


food 


H#015 


19 55 26.8 


-55 09 20 


3 


gl955268-550920 


0.045 


16.49 


1.43 


4.30 


3.4 


good 


H#005 


19 56 51.4 


-54 58 37 


1 


gl956514-545837 


0.061 


16.65 


1.37 


4.35 


5.2 


bad 


H#010 


19 56 51.6 


-5519 56 


1 


gl956516-551956 


0.056 


16.35 


1.23 


4.23 


5.3 


good 


H#004 


19 56 56.7 


-55 47 30 


1 


gl956567-554730 


0.018 


14.67 


1.03 


3.62 


8.1 


good 


H#014 


195711.2 


-55 25 09 


1 


gl957112-552509 


0.059 


16.40 


1.33 


4.41 


4.4 


good 


H#009 


1957 22.2 


-55 0814 


1 


gl957222-550814 


0.016 


15.36 


1.25 


4.31 


10.6 


good 


H#011 


195733.7 


-55 34 41 


1 


gl957337-553441 


0.017 


14.45 


1.26 


4.32 


8.1 


good 


H#001 


19 58 00.3 


-55 33 29 


1 


gl958003-553329 


0.056 


16.90 


1.46 


4.55 


3.4 


good 


H#012 


19 5812.8 


-55 40 53 


9 


gl958128-554053 


0.017 


15.62 


0.16 


4.07 


14.5 


good 


H#006 


19 58 29.1 


-55 0913 


1 


gl958291-550913 


0.058 


16.14 


1.40 


4.65 


4.9 


bad 


H#013 


19 58 45.0 


-55 3511 


1 


gl958450-553511 


0.058 


16.11 


1.24 


4.20 


5.5 


good 


H#002 


19 5617.6 


-55 0711 


1 


gl956176-550711 


0.060 


16.96 


1.73 


3.97 


2.4 


bad 



1. J-band effective radii from 2MASS XSC. 

2. The 'bad' flag indicates the three hexabundles that were damaged during the run. 



under-estimate the true flux in the star by as much as 10 
percent. This can be avoided by performing a more careful 
PSF extraction of the star over the relevant FoV, but in this 
case we used the simpler method and this should be borne 
in mind for the following results. The summed spectrum was 
then compared to the tabulated values for the spectropho- 
tometric standard. 

The results of this analysis are shown in Figures [9] 
and [To] for the blue and red arms of the spectrograph respec- 
tively. These throughput curves have not been corrected for 
telluric features, which are still visible. The curves therefore 
show the measured throughput of the entire system from at- 
mosphere to detector, including the AAT and the AAOmega 
spectrograph. Of note is the fact that the throughput falls off 
quite quickly in the blue arm. This is not ideal but the source 
of the effect is known. The fibre train for SAMI consists of 
~42m of AFS105/125Y fibre matched to the fibre used in 
the maimfacture of the hexabundles themselves. This type of 
fibre has poorer throughput in the bl ue than the Polyrn icro 
FBP fibre used, for example, in 2dF (|Smith et al.ll2004 ). 

Throughput values measured from the standard star 
were compared with the lab-tested throughput from Sec- 
tion [331 This is a difficult comparison because the lab tests 
use Bessel _B-band and _R-band filters and the through- 
put measured from the standard star observations changes 
rapidly through the bandpass of the Bessel _B-band filter. 
The lab tests were corrected for the known throughput in the 
B- and i?-bands of the atmosphere (0.72, 0.89 for B and R 
respectively), telescope including corrector (0.77, 0.80), and 
spectrograph including CCD (0.17, 0.33). The lab tests have 
uncertainties due to fibre-to-fibre differences. The data re- 
duction includes fibre-to-fibre throughput corrections for the 
standard stars, which scales the throughput of each individ- 
ual fibre to be closer to that of the fibre in the centre of a rib- 
bon. We therefore compare the throughput measured from 
the observations to core 6 in Table [1] The SAMI standard 
star images had a steep variation in observed throughput 
across the _B-band ranging from 26-46 percent (for SAMI 
only), while the lab-tests found 32 ± 7 percent. The 7?-band 



lab measurement is not directly comparable to the on-sky 
measurement, as the _R-band filter used in the lab has a peak 
transmission wavelength outside the range of the AAOmega 
lOOOR grating. However, we show an indicative comparison 
between direct and lab-based throughput measurements in 
Fig. [To] The throughput comparisons, while limited by the 
differences in beam shape and filters used, shows that there 
are no other significant losses in the fully installed system 
that are not accounted for in the above analysis of through- 
put. 

From the measured throughput of the SAMI system 
we can calculate the limiting surface brightness for typical 
observations, assuming photon counting errors. Here we as- 
sume 3 hours exposure time per field and a dark sky {V and 
R band sky brightnesses of 21.5 and 20.8 mags arcsec"^). 
For a S/N of 5 the estimated limiting surface bright- 
ness is 22.9 and 22.5 magarcsec"'^ for the V and R bands 
respectively (on the Vega scale). For a fiducial galaxy sur- 
vey sample, which is flux limited to r < 16.5 and redshift 
z < 0.13, these surface brightnesses limits will allow the cur- 
rent SAMI system to achieve S/N> 5 at 1 re for 80 percent 
of galaxies. 

5.4 Galaxy observations 

One galaxy held with thirteen galaxies was observed dur- 
ing the commissioning run. Three 40-minute exposures were 
taken at this pointing on the night of 2 July. These were re- 
duced using 2dFdr and stacked, resulting in a total of 2 hours 
integration on this field. A spectrophotometric standard star 
was observed on the same night as the objects and was re- 
duced in an identical way. The extracted star spectrum was 
then used to correct the galaxy data for instrument through- 
put, though the data were not fully fiux-calibrated. The de- 
tails of the 13 target galaxies are given in Table [3] 

Figure [TT1 shows the ten galaxies from the observed field 
that had usable data. The remaining 3 hexabundles were 
damaged during commissioning. The images were created 
by summing over a wide wavelength range (6300-7300 A in 
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Figure 11. IFU images of the observed SAMI galaxies obtained by summing over the spectrum in each spaxel. The X and Y positions 
are in arbitrary pixel coordinates. Each spaxel is represented by a circle with a 10-pixel diameter corresponding to 1.8". The colour scales 
indicate the logarithm of the flux in arbitrary units. 
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Figure 12. Red spectra for the galaxy observed writh bundle 
H#009. The top plot shows the integrated spectrum converted to 
rest frame co-ordinates. The central plot shows the spectrum from 
core 1 of the hexabundle and the bottom plot shows the spectrum 
from core 28. The Ha, [N II]AA6548, 6583 A and [S II]AA6716, 
6731 A emission lines are clearly visible. 



the observed frame); the colour bars indicate the logarithm 
of the flux in arbitrary units. The discs and bulges of the 
larger galaxies are readily apparent 



Figure 13. Blue spectra for the galaxy observed with bundle 
H#009. The top plot shows the integrated spectrum converted 
to rest frame co-ordinates. The central plot shows the spectrum 
from core 1 of the hexabundle and the bottom plot shows the 
spectrum from core 28. The H/3and [O III]A5007 A are clearly 
visible in the integrated and central spectra but have lower S/N 
in the outer spectrum. The masked regions in the spectra are due 
to bad pixels in the AAOmega blue CCD. 
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Figure 14. A BPT diagram showing the positions of the three 
SAMI galaxies for which it was possible to fit all lines. Also 
shown are the theoretical maximum for star-forming galaxies 
(blue curve) and the empirical cut-off for star-forming galaxies 
(red curve); AGN and LINERs are separated, empirically, by the 
black line (AGN above, LINERs below); see text for details. The 
green dots show the positions of each of the individual spaxels for 
the galaxy observed with bundle H#009. 



Representative spectra of the galaxy observed in bundle 
H7^009 are shown in Figures [12] and 1131 The three spectra 
in each figure correspond to the integrated spectrum of the 
galaxy (top panel), the spectrum observed in the central core 
(core 1; the central panel) and the spectrum observed in an 
outer core (core 28; the bottom panel). The integrated spec- 
tra are simply the sum of all the individual spaxel spectra in 
a particular IFU. Various emission and absorption features 
are seen in the galaxy spectra. In the red wavelength range 
the Ha, [N II]AA6548, 6583 A and [S II]AA6716, 6731 A emis- 
sion lines are clearly visible in all three spectra, though with 
lower S/N in core 28. In the blue wavelength range we see 
the and [O III]A5007 A emission lines, with an H/3 ab- 
sorption feature also visible. The data in core 28 have lower 
S/N than the central core, but still sufficient for analysis. 

In examining the results, we first studied the integrated 
spectrum of each galaxy. Seven galaxies show medium to 
strong Hq and [N II] emission (along with the [S II] doublet 
at 6716 and 6731 A) and three of these also show the H/3 
and [O III] lines. 

The ratios of these lines are often used to clas- 
sify objects by means of a so- called BPT diagram 
jBaldwin. Phillips fc TerlevichI [l98ll 'l . Each line was fitted 
with a Gaussian to measure the line strength. In the case of 

a weak absorption trough was observed in the summed 
spectra of all three galaxies; this trough was fitted separately 
and accounted for when calculating the H/3 line strength. 
Figure [14] shows the galaxies on a BPT diagram (large di- 
amonds). In Figure [14] the blue curve shows t he theoretical 
maxim um line for star-forming galaxies after iKewlev et al.l 
while the red curve shows the empiric al cut-off for 
star-forming galaxies from 'Kauffmann et al." ('200 j). The 
black line is an empirical line, also from Kauffma nn et all 
(|2003l ). that separates AGN (above the line) from LINERs 
(below the line). 

The galaxy H#009 (an image of which is shown in the 



left-hand panel of Figure fTSj) was then studied in more detail, 
using the full spatially-resolved set of spectra. The line ratios 
discussed above were calculated for each individual spaxel 
in the IFU FoV, yielding spatially-resolved line ratio maps, 
which are shown in the centre and right panels of Figure [151 
There is a systematic spatial trend in both line ratios moving 
away from the plane of the galaxy's disc. The green circles 
in Figure [14] shows the line ratios from each spaxel plotted 
on a BPT diagram. A strong trend is seen in the sense that 
the spaxels further from the plane of the disc have a harder 
ionising mechanism, refiecting the spatial trend seen in the 
line ratio maps. A single-fibre spectrum of this galaxy would 
have seen only the central few arcseconds of the object and 
resulted in classification as a normal star-forming galaxy. It 
is clear from the SAMI results, however, that there is a more 
interesting story behind this galaxy, which will be discussed 
in detail in a future paper (Fogarty et al. in preparation). 

Gas kinematics were examined for six of the galaxies us- 
ing the Ha emission line. In each spaxel of each galaxy the 
Ha line was fitted by a Gaussian function and the kinematic 
properties (mean velocity and velocity dispersion) extracted 
from the fit. The resulting velocity maps are shown in Fig- 
ure [16] The colour bars show velocity in kms~^. From these 
maps it is clear the galaxies H#009, H#010, H#012 and 
H#016 are fairly typical rotating discs; the velocity maps 
for H#001 and H#004 are possibly more complex. 

From the wealth of spatially-resolved information pre- 
sented for the handful of galaxies studied in this single field 
during the commissioning run, it is clear that SAMI can ex- 
plore a very wide parameter space with enormous science 
potential. 



6 CONCLUSIONS 

In this paper we have presented a new instrument, the 
Sydney- AAO Multi-object IFS (SAMI). SAMI makes use of 
astrophotonic technology in the form of hexabundles (multi- 
core fibre bundles) to enable simultaneous IFU observations 
of 13 objects over a 1-degree diameter accessible field. Each 
IFU contains 61 elements, each 1.6 arcsec in diameter, giving 
a FoV for each IFU of 15 arcsec. SAMI has now been com- 
missioned on the Anglo- Australian Telescope (AAT) and we 
demonstrate its science potential via preliminary observa- 
tions of galaxies selected from the 6dF Galaxy Survey. 

We make the case that multiplexed integral field spec- 
troscopy is the natural next step in galaxy surveys, which 
to date have been dominated by multiplexed single-aperture 
observations. The extra information gained by IFU spec- 
troscopy, combined with the statistical power of a large sur- 
vey, will enable a fundamental step forward in our under- 
standing of galaxy formation and evolution by distinguishing 
the spectroscopic properties of the major structural compo- 
nents of galaxies. Assuming that a survey could observe 3 
fields per night (i.e. nominal exposure times between 2 and 
3 hours), SAMI will allow 10,000 galaxies to be targeted in 
260 clear nights on the AAT. 

Much greater gains could be made with an increased 
FoV and larger numbers of IFUs. A key requirement for 
making full use of the larger fields of view available on cur- 
rent and future telescopes is to have sufficient spectrographs 
to handle all the fibres from the IFUs. A 2-degree diame- 
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Figure 15. Analysis of the spatially-resolved spectra of the galaxy observed with H#009. The left-most pa nel shows the field o f the 
hexabundle (red circle, diameter 15 arcsec) superimposed on a UKST bj-band image from SuperCOSMOS l|Hamblv et al.ll200lf l: the 
image is 1 X 1 arcminute in size. The other two panels show the spatially- resolved line ratio maps for the galaxy, with [N II]A6583/Ha 
in the centre and [O III] A5007/H/3 to the right. 
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Figure 16. Gas kinematics for six of the SAMI galaxies. On the left is a SuperCOSMOS jHamblv et al.ll200 bj-band image of each 
galaxy, with the IFU FoV represented by the red circle. Each image is ~1 arcmin on a side. The SAMI Ha velocity maps arc shown on 
the right hand side, with the X and Y coordinates in pixel units. 



ter FoV (e.g. that available with the AAT's 2dF corrector) 
could accommodate ~50 or more IFU bundles comparable 
to those in SAMI; this would, however, require ^^3300 fibres 
(including sky fibres) . Such large numbers of fibres naturally 
drive designs towards mass-produced fixed-format spectro- 
graphs, such as those currently being developed for projects 



such as MUSE (|Bacon et al.ll2004l ) and VIRUS (|Hill et all 
120061 ). 

More development work on the bundle technology is 
under way. One possible improvement is to force the fi- 
bre cores into a perfect hexagonal configuration to sim- 
plify the data analysis, particularly when using dithering 
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to smooth out the discrete samphng. This will not be easy 
to do with our current core size while at the same time en- 
suring an AR-coated polished front facet that is not overly 
stressed by the hexagonal grid. With the successful commis- 
sioni ng on the AAT of fibre-based OH suppress i on technol- 
ogy (iBland-Hawthorn. Englund fc Edvellf l20oi : lEUis et all 
|2010| ). we can envisage whole bundles that are fully sky- 
suppressed at near-infrared wavelengths. Once the near- 
infrared sky is rendered as dark as the optical sky, exposure 
times in both regimes are expected to be the same and it 
makes sense to run optical and infrared spectrographs si- 
multaneously. 
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